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Preface 


This  thesis  Is  part  of  a continuing  study  of  the 

problem  of  computer  analysis  of  alpha  particle  spectra. 

The  performance  of  the  program  under  Ideal  conditions  was 

210 

studied  using  Gaussian  curves.  Po  was  used  to  prepare 
sources  of  different  thicknesses  to  obtain  a library  of 
curves.  These  were  then  used  to  analyze  samples  of  americium 
and  plutonium.  Using  a reference  peak  with  a full  width  at 
half  maximum  that  matched  the  sample  gave  the  most  accurate 
analysis. 

Several  people  contributed  to  the  success  of  this 
thesis.  Dr.  George  John  and  Dr.  Richard  Hagee,  my  thesis 
advisors,  provided  ideas,  resources,  and  encouragement  all 
summer.  Harold  Kirby  and  Roland  Armani  contributed  their 
ideas  and  experience  in  the  preparation  of  sources. 

Dr.  Dave  Hardin  and  Captain  Bruce  Stinson  were  both  helpful 
in  my  programming  efforts.  To  them  in  particular,  and  to 
all  of  the  AFIT  physics  laboratory  and  metal  shop  personnel 
as  well,  I extend  my  sincere  thanks. 

Finally,  I would  like  to  thank  Mrs.  Donna  Hadley,  whose 
typing  contributed  significantly  to  the  successful  completion 
of  this  thesis. 
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Abstract 

Sources  of  Po^®  of  varying  thicknesses  were  prepared 
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by  an  evaporative  precipitation  method.  Sources  of  Am 
240 

and  Pu  were  analyzed  with  the  spectra  of  these  sources 
as  reference  peak  functions  using  a modification  of  an 
existing  computer  code.  The  results  showed  a reduction  of 
residual  from  using  a reference  peak  with  a full  width 
at  half  maximum  close  to  that  of  the  analyzed  source.  The 
performance  of  the  computer  code  was  also  studied  using 
randomized  Gaussian  peaks. 


COMPUTER  CODE  TO  ANALYZE  ALPHA 
SPECTRA  USING  A SPECTRAL 
STRIPPING  APPROACH 

I . Introduction 

This  report  is  part  of  a continuing  study  of  alpha 
spectra  obtained  from  solid  state  detectors.  Ii*  general 
purpose  of  this  study  is  to  learn  how  to  use  these  spectra 
to  obtain  a more  accurate  analysis  of  small  quantities  of 
alpha-emitting  nuclides.  The  specific  purpose  of  this 
thesis  was  to  take  a computer  program  which  analyzes  alpha 
spectra  using  a spectral  stripping  approach  and  modify  it  to 
compensate  for  the  differences  in  the  shape  of  spectra 
obtained  from  different  sources. 

Basic  Problem 

Solid  state  alpha  spectroscopy  is  a technique  widely 
used  to  analyze  microscopic  quantities  of  unknown  substances 
to  determine  the  amounts  in  the  sample  of  various  heavy 
nuclides  such  as  plutonium.  However,  the  accuracy  of  this 
analysis  is  limited  by  several  problems  inherent  in  alpha 
decay  and  its  detection.  The  spectrum  from  particles  of  a 
single  energy  is  spread  over  the  range  of  energies  from  the 
initial  energy  of  decay  practically  down  to  ze^o.  This  can 
be  seen  clearly  in  Figure  1,  which  is  a full-range  spectrum 
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from  a source  of  Radium  DEF  (Po  ).  In  addition,  most 
alpha  emitters  can  decay  to  more  than  one  energy  state  In 
the  daughter  nucleus  yielding  alpha  particles  of  different 
discrete  energies  in  tho  process.  The  spectra  from  the 
different  energies  are  then  overlapped,  as  seen  in  Figure  2, 
which  represents  part  of  a spectrum  from  Americium  241. 
Finally,  the  shape  of  the  curves  from  the  individual  alpha 
energy  peaks  and  of  the  resulting  combined  spectra  varies 
from  sample  to  sample  because  of  differences  in  the  self- 
absorption.  This  is  shown  in  Figures  3 and  4,  which  are 
overlapping  normalized  spectra  from  different  samples  of 
Po^*®  (Figure  3)  and  Pu^®  (Figure  4). 

If  a sample  is  to  be  analyzed  accurately,  the  counts  of 
a given  energy  must  be  separated  into  the  individual  com- 
ponents from  the  particles  of  different  original  energies. 
This  requires  a reference  curve  of  the  correct  shape  for 
the  particular  sample  to  be  analyzed.  The  components  of 
the  entire  energy  range  for  each  energy  alpha  particle  are 
then  summed.  Finally,  the  contributions  from  all  the  alpha 
particles  for  each  nuclide  are  summed  to  give  the  total 
number  of  decays  for  that  particular  nuclide. 

Background 

This  study  was  specifically  designed  to  be  a continu- 
ation of  work  done  by  R.  J.  Hartley  (Ref  5),  He  had  prepared 
Program  ALPHAFIT  with  the  help  of  Dr.  Philip  Poirer  from 
Program  AUTOFIT  which  was  written  by  J.  R.  Comfort  of 


FIG  3 . Rq  OEF  * 1»  ( UNOXIDIZcD ) 0.6  KEV/CHRNNcL  40.000  SEr 
RR  OEF  * 1J  (UN0XID’2ED)  O.C  KcV'PHRM^EL  40s000  ScT 


FIG  4 . PU  240  * 1 flSMflSl  77  0.6  KcV/CH4NNEL  20.000  SEC 

PU  240  * 2 77  0.6  KcV/CHflNNEL  10.000  SEf 


Argonne  National  Laboratory.  AUTOFIX  had  been  written  to 

analyze  gamma  ray  spectra  and  Hartley  modified  it  to  analyze 

alpha  particle  spectra.  He  then  studied  its  behavior  using 

computer  generated  Gaussian  peaks  and  noted  an  anomaly  in 

the  computed  error,  which  increased,  dropped,  and  then 

Increased  again  as  the  two  peaks  in  the  spectrum  to  be 

analyzed  were  set  progressively  closer  together.  He  also 

24 1 240 

used  it  to  analyze  samples  of  Am  and  Pu  using  the 
210 

spectrum  from  Po  as  a reference  curve  shape.  However, 
the  results  obtained  showed  large  fluctuations  in  both  the 
computed  errors  and  the  computed  peak  heights  when  either 
the  sample  of  the  same  analyzed  nuclide  or  the  sample  used 
as  a reference  was  changed.  These  fluctuations  were 
believed  to  be  related  to  the  variations  in  the  source 
thicknesses. 

Problem  Statement 

The  problem  is  to  develop  a computer  program  that  will 
take  a spectrum  from  a multichannel  analyzer,  separate  it 
into  its  separate  components,  and  compute  the  quantities 
of  the  various  nuclides  present. 

Plan  of  Attack 

The  basic  approach  selected  was  to  modify  Program 
ALPHAFIT  to  use  a library  of  curves  with  shapes  represen- 
tative of  a range  of  source  thicknesses.  The  program  would 
then  choose  the  curve  with  characteristics  matching  those 
of  the  source  to  be  analyzed  and  use  that  curve  as  a 


reference  to  analyze  the  spectrum  of  the  unknown  sample. 

A set  of  reference  samples  were  to  be  prepared  using  either 
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Radium  DEF  or  just  Polonium  and  reference  spectra 
obtained  from  these  samples.  Also,  the  anomaly  in  the 
computed  errors  was  to  be  studied  and  the  program  modified 
if  necessary  to  correct  the  problem. 


Summary  of  Work  Done 

The  program  was  first  modified  to  read  data  directly 

from  magnetic  tape  or  permanent  disc  storage.  In  order  to 

reduce  the  core  memory  required  and  shorten  the  turnaround 

time,  the  program  was  changed  from  analyzing  20  peaks  of 

2048  channels  each  to  analyzing  5 peaks  of  512  channels. 

Subroutines  to  generate  Gaussian  peaks  were  then  written  and 

the  program's  behavior  studied. 

At  the  same  time,  a series  of  samples  of  Ra  DCF  was 

prepared  with  varying  quantities  of  lead  nitrate  added  to 

give  samples  of  the  same  activity  but  varying  thickness.  An 

evaporative/precipitation  technique  developed  by  H.  W.  Kirby 

of  the  Mound  Facility  in  Miamisburg,  Ohio,  was  used  to 

prepare  these  samples.  An  electroplating  technique  was  also 

tried  but  sufficient  time  was  not  available  to  perfect  it. 

The  spectra  from  these  samples  were  then  recorded  on  magnetic 

tape  and  used  to  analyze  spectra  from  existing  samples  of 
241  240 

Am  and  Pu  . However,  the  range  of  thicknesses  was  only 
appropriate  for  one  of  the  plutonium  samples.  This  sample 


was  analyzed  using  the  reduced  version  of  the  program  and 


9 


analyzed  again  after  the  program  was  converted  back  to 
analyzing  2048  channels  for  up  to  10  peaks.  The  results  of 
these  analyses  confirmed  the  value  of  the  library  of 
curves  approach. 


Overview  of  Thesis 

In  this  thesis,  the  theory  pertinent  to  the  analysis 
of  alpha  spectra  is  first  discussed.  Ihen,  the  equipment 
and  procedures  used  are  described.  The  results  of  the 
experiments  are  presented  and  analyzed.  The  changes  to  the 
program  are  explained.  The  results  of  the  computer  analyses 
are  presented  and  discussed.  Finally,  the  conclusions 
drawn  from  this  study  are  set  forth  along  with  recommendations 
for  future  work. 


t 


II.  Theory 


l! 

The  analysis  of  Che  amount  of  different  nuclides  by  use 
of  alpha  speccroaecry  requires  an  understanding  of  the 
nature  of  the  alpha  spectra  emitted  by  different  nuclides, 
the  processes  which  perturb  the  energy  of  the  alpha  particles 
emitted , and  the  response  of  the  detector  to  the  alpha 
particles.  These  are  discussed  in  this  chapter. 

Alpha  Decay 

Alpha  decay  occurs  as  a transition  from  the  nucleus  of 
one  of  the  heavy  Isotopes  to  a lower  energy  state  in 
another  isotope  by  the  expulsion  of  a helium  nucleus.  Each 
parent  nuclide  emits  alpha  particles  of  certain  specific 
energies,  as  listed  in  Table  I.  The  exact  energy  released 
is  determined  chiefly  by  the  quantum  state  of  the  resulting 
daughter  nucleus.  Once  emitted,  the  alpha  particle  may  lose 
varying  amounts  of  energy  by  a number  of  processes  up  to  and 
including  detection. 

If  the  emitting  source  material  is  extremely  thin,  some 
alpha  particles  may  leave  the  source  without  losing  any 
energy  in  collisions.  However,  if  the  source  is  even  a few 
atoms  thick,  most  alpha  particles  will  lose  some  energy  in 
ionizing  collisions  with  the  atomic  electrons  of  the  source 
material.  The  energy  lost  Increases  with  the  path  length 
through  the  material  which  is  the  thickness  of  the  material 

10 


I so tope 

Type 

Decay 

Half-Life 

Percent 
of  Decays 

Pb210(RaD) 

3* 

20.4  years 

0.061  (max) 

100 

Bl210(RaE) 

3" 

5.013  days 

1.160  (max) 

100 

Po210(RaF) 

a 

138.40  days 

5.305 

100 

*»241 

a 

458  years 

5.486 

86 

5.443 

12.7 

5.389 

1.3 

Pu240 

a 

6580  years 

5.168 

76 

5.123 

24 

Pu239 

a 

2.44 

x 10^  years 

5.157 

73.3 

5.145 

15.1 

5.107 

11.5 

divided  by  the  sine  of  the  angle  between  the  path  and  the 
plane  of  the  source  surface.  This  energy  loss  Is  termed 
self -absorption. 

A few  alpha  particles  will  Initially  be  emitted  away 
from  the  detector  but  will  collide  with  a nucleus  in  the 
source  Itself  or  the  source  backing  material.  After  a 
large  change  in  direction  and  a large  loss  of  energy  in  the 
collision,  some  of  these  particles  will  reach  the  detector 
and  deposit  their  remaining  energy.  This  is  termed  back- 
scatter. 

Any  material  such  as  air  between  the  source  and  the 
detector  will  absorb  energy  from  the  alpha  particles,  so 
alpha  spectrometry  should  be  conducted  in  an  evacuated 
chamber. 

When  they  reach  the  detector,  the  alpha  particles  may 
have  to  pass  through  a dead  layer  at  the  surface  of  the 
detector  in  which  the  electrons  freed  by  collisions  do  not 
get  collected  to  contribute  to  the  pulse  of  current  from 
the  detector.  As  in  self -absorption,  the  energy  lost  is 
dependent  on  the  thickness  of  the  dead  layer  and  the  angle 
of  the  particle's  path  through  it. 

Finally,  the  alpha  particle  reaches  the  active  region 
of  the  detector,  where  it  will  lose  its  energy  in  collisions 
with  electrons.  These  electrons  produce  a pulse  of  current 
which  is  amplified,  shaped,  and  counted  in  a multi-channel 
analyzer. 


Alpha  Particle  Interaction 

The  Interaction  of  alpha  particles  with  matter  is 

strongly  influenced  by  two  characteristics  of  the  alpha 

particle i its  enormous  mass  relative  to  an  electron  mass 

(M  « 7000  iO  and  its  double  positive  charge.  The  mass 
ct  e 

causes  the  alpha  particle  to  travel  in  a virtually  straight 
line  except  for  Rutherford  scattering  with  another  nucleus. 
It  also  gives  alpha  particles  a relatively  low  velocity  for 
their  high  energy;  for  example,  an  alpha  particle  of  energy 
5 MeV  has  a speed  only  37.  that  of  light.  An  electron  of 
the  same  energy  would  travel  at  99%  the  speed  of  light. 

Thus  an  alpha  particle  has  thirty  (or  more)  times  as  long 
as  an  electron  for  its  double  charge  to  interact  with  an 
electron  from  each  atom  it  passes.  Alpha  particles  are 
easily  stopped  by  only  a thin  layer  of  matter  and  may  lose 
considerable  energy  before  they  even  leave  their  source 
material . 

The  stopping  power  or  specific  energy  loss  dE/dX  is 
given  approximately  by  the  following  equation  from  Evans 
(Ref  40i637)i 


dX 


4,eV  \ ° 

—pr  GAW 


In 


2m0V2 

I 


fn(  1 - e2) 


(1) 


with 

e ■ electronic  charge 
m^  * electron  rest  mass 
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E - kinetic  energy  of  the  primary  (alpha)  particle 
ze  » charge  of  the  primary  particle 
V * velocity  of  the  primary  particle 
■ Avogadro ' s number 

p ■ density  of  the  material  in  grams  per  unit  volume 

GAW  = the  gram  atomic  weight  of  the  element 

Z ■ atomic  number  of  the  absorber 

8 “ V/C,  where  C = velocity  of  light  in  a vacuum 

I = geometric -mean  ionization  and  excitation  potential 
for  the  absorbing  atoms.  This  cannot  be  calculated 
accurately  and  is  usually  considered  to  be  constant 
for  each  element.  It  must  be  determined  experi- 
mentally for  each  element. 

It  is  frequently  more  convenient  to  determine  the  loss 
per  differential  mass/area  thickness  traveled  through.  This 
is  done  by  dividing  both  sides  by  p,  the  density.  For  a 
substance  with  several  different  elements  such  as  lead 
nitrate,  PbCNO^^,  it  is  necessary  to  compute  the  energy  loss 
for  each  element  separately  and  then  add  the  individual 
losses. 

Alpha  Particle  Detection 

This  study  deals  with  the  detection  of  alpha  particles 
by  solid  state  surface  barrier  detectors.  When  an  alpha 
particle  enters  the  active  region  of  the  detector,  the 
electrons  with  which  it  interacts  and  the  holes  left  behind 
them  move  under  the  influence  of  the  applied  bias  voltage. 

This  gives  a transient  current  which  can  be  converted  into  a 


voltage  pulse  and  counted  by  the  associated  electronic 
equipment.  The  height  of  this  pulse  can  be  expressed  by  the 
following  equationi 

ph  * [(Edep/w)  x e x y]/C  = Nev/C  = q/C  (2) 

where 

Edep  = energy  deposited  by  the  alpha  particle 

w * average  energy  required  for  the  production  of 
an  electron-hole  pair  (eV/e-h  pair) 

e ■ electron  charge  (coulombs) 

Y ■ charge  collection  efficiency.  This  is  essen- 
tially . 1 in  detectors  made  of  high  quality 
single  crystals. 

N **  Edep/W  - number  of  electron-hole  pairs 
produced 

q = N€y  = total  charge  collected  (coulombs) 

C ■ capacitance  of  the  detector  (farads) 
ph  **  pulse  height  in  volts  at  the  detector 

These  pulses  are  then  amplified  and  shaped  by  the 
preamplifier,  linear  amplifier  and  biased  amplifier  and 
transmitted  to  a multichannel  analyzer  which  counts  the 
number  of  pulses  with  height  between  h and  h + Ah  during 
the  acquisition  time.  These  are  collected  in  a number  of 
bins  to  form  a pulse-height  distribution. 

Because  of  the  statistical  nature  of  the  energy  loss 
process,  the  pulses  will  be  distributed  over  a range  of 
channels  or  energies  even  if  the  alpha  particles  have 


precisely  the  same  energy  when  incident  on  the  detector 
face.  This  statistical  spreading  is  not  caused  by  a distri- 
bution of  either  the  original  alpha  energy  or  the  energy 
deposited  in  the  detector.  It  is  caused  by  statistical 
variations  in  the  number  of  electron-hole  pairs  created  by 
each  alpha  particle  as  it  deposits  its  energy  in  the 
detector.  Additional  spread  in  the  pulse  height  distribution 
is  caused  by  a spread  in  the  energy  of  the  alpha  particles 
caused  by  self -absorption  in  the  source  material,  the  dead 
layer  at  the  surface  of  the  detector,  and  variations  in 
the  angle  at  which  the  alpha  particle  passes  through  the 
dead  layer  and  the  source.  Further,  though  minor,  contri- 
butions to  the  total  arise  from  large  angle  scatters  from 
the  backplate  of  the  source  and  the  surroundings.  These 
combine  to  give  the  tailing  effect  on  the  low  energy  side 
of  the  peak. 

Electronic  noise  in  the  electronics  also  gives  a 
slight  increase  in  the  spread  but  this  will  be  the  same  for 
both  the  reference  source  and  the  unknown  when  they  are 
analyzed  with  the  same  detector  system. 

Relation  of  Counts  to  Activity 

Each  unknown  spectrum  is  analyzed  to  determine  the 
number  of  peaks  and  the  area  of  each  peak.  This  area  is 
the  number  of  counts  per  channel  summed  over  all  of  the 
channels.  This  is  related  to  the  source  activity  by  the 
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following  formula i 


Nt/tc  - A x G x £ss  x fK  x E x £t  x fa 


where 

Nt  - total  number  of  counts  = the  combined  areas 
for  all  the  peaks  of  each  isotope 

tc  ■ the  counting  time 

A - the  source  activity 

G ■ solid  angle  subtended  by  the  active  layer  of 
the  detector  relative  to  the  source 

fao  ■ surrounding  scatter  factor  (negligible  for 
s s 

a particles) 

f ■ window  factor,  that  is,  absorption  losses 
(generally  * 1 for  zero  losses) 

E - intrinsic  efficiency  of  the  detector  in  counts/ 
alpha  (generally  1) 

ft  «=  dead  time  correction  (generally  1 for  zero 
dead  time) 

f ■ number  of  alphas  per  nuclear  transformation 


Need  for  Computer  Peak  Analysis 

In  alpha  spectrometry,  the  energies  and  relative  num- 
bers of  alpha  particles  emitted  by  each  nuclide  are  assumed 
to  be  known  from  other  studies,  the  results  of  which  are 
compiled  and  tabulated  in  works  such  as  Lederer  (Ref  7). 
These  works  provide  an  outer  bound  on  the  set  of  possible 
nuclides,  which  can  frequently  be  further  restricted  by 
other  information  about  the  unknown  source.  If  only  one 


nuclide  is  present,  then  all  of  the  nonbackground  counts  are 
from  that  nuclide  and  computer  analysis  is  unnecessary.  It 
is  only  in  those  cases  where  the  relative  numbers  of 
nuclides  are  unknown,  while  the  possible  nuclides  and 
energies  and  proportions  of  alphas  emitted  by  each  of  those 
nuclides  are  known,  that  computer  peak  analysis  is  both 
necessary  and  feasible. 

Computer  Peak  Analysis 

Each  spectrum  is  represented  by  an  array  [A^]  which  is 
the  sum  of  a set  of  N components  Z^(xn),  one  for  each  of  the 
N peaks.  Each  element  in  A is  the  sum  of  the  corresponding 
elements  of  each  component 

Ai  “ = . an  Zi(*n>  <4> 

n*i 

where  Z^(xR)  represents  the  counts  in  channel  i from  a 
source  component  that  has  its  peak  located  at  xn>  Ihe 
locations,  xn,  and  relative  strengths,  an,  are  varied  until 
a satisfactory  fit  of  the  test  array  £a]  to  the  data 
array  [y]. 

Mathematically,  this  is  accomplished  by  forming  a 
chi -squared  function 

N 7 7 

F - £ (Yt  - A.r/CUT.r  (5) 

i-1  1 1 1 

which  is  summed  over  the  region  being  analyzed.  The 
weighting  function  UY^  is  + Yt  , where  Y^  is  the 


ordinate  of  the  ith  data  point,  that  is,  the  counts  in 
channel  i.  The  inclusion  cf  10  slightly  increases  the 
weighting  of  the  peak  region  relative  to  the  tail  and  back- 
ground and  avoids  division  by  zero. 

Since  Poisson  statistics  adequately  describe  the 

statistical  fluctuations  of  the  number  of  counts  in  each 

2 

channel,  the  weighting  function  (U Y^)  is  approximately 
equal  to  the  variance.  Thus,  the  chi- squared  function 
represents  the  ratio  of  the  observed  spread  in  the  analyzed 
data  to  the  expected  spread  or  variance. 

The  xn's  are  determined  by  computing  the  function  for 
locations  to  either  side  of  xn  and  shifting  the  xn's  in 
the  direction  that  reduces  F.  When  a local  minimum  is 
found  or  when  directed  by  the  input  data,  that  particular 
%n  is  fixed. 

The  an's  are  then  determined  from  direct  solution  of 
Eq  (4)  by  matrix  inversion  techniques. 
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Ill . Experimental  Equipment  and  Procedures 


Introduction 

This  section  covers  the  techniques  and  equipment  used 
for  alpha  spectroscopy  and  the  preparation  and  assay  of  the 
reference  sources  of  Radium  DEF. 

Alpha  Spectroscopy  Equipment 

The  equipnent  used  in  these  experiments  consists  of  a 
vacuum  chamber,  surface  barrier  detector,  power  supplies 
and  amplifiers,  and  multi-channel  analyzers. 

Vacuum  Chamber.  The  vacuum  chamber  was  a steel 
cylinder,  10  cm  in  diameter  and  15  cm  in  height  (See 
Figure  5).  The  cylinder  is  open  at  the  top  and  is  mounted 
around  an  0-ring  clamp  to  the  top  plate.  There  are  four 
openings  through  the  top  plate.  Two  of  these  lead  to  gas 
valves.  One  valve  is  connected  to  the  vacuum  pump,  while 
the  other  can  be  opened  to  the  room  air.  The  third  is  the 
mount  and  connection  for  the  detector.  A BNC  connector 
above  the  plate  is  electrically  connected  through  the  plate 
to  a threaded  mount  or  microdot  connector  into  which  the 
detector  is  screwed.  A BNC  ”L"  connector  leads  directly 
from  this  to  the  preamplifier  which  is  supported  by  a clamp 
and  stand.  The  fourth  opening  is  an  O-ring  fitting  for  a 
thin  rod  which  is  attached  to  a 3 cm  x 3 cm  horizontal  plate 
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Fig.  5.  Raised  and  Parallel  Views  of  Vacuum  Chamber  for 
Alpha  Spectroscopy  (Not  to  Scale) 


2! 


f 


to  hold  samples.  The  distance  from  source  to  detector  Is 
varied  by  sliding  this  rod  up  or  down. 

2 

Detector.  An  ORTEC  100  mm  surface  barrier  detector 
(SN  8-907A)  with  a raicrodot  mount  was  used.  This  detector 
takes  a maximum  bias  of  50.0  volts.  A gross  bias  of 
50  volts  was  applied  for  a net  bias  of  30-35  volts. 

Amplifiers.  The  output  from  the  surface  barrier 
detector  was  fed,  in  turn,  to  an  ORTEC  Model  109-A  preampli- 
fier, a Tennelec  Model  TC203BLR  linear  amplifier,  and  an 
ORTEC  Model  408  biased  amplifier. 

Multi-channel  Analyzers.  Two  multi-channel  analyzers 
were  used  to  convert  the  output  pulses  into  a count  spectrum. 
An  ND100  was  used  to  obtain  512  channel  spectra  for  prelim- 
inary observations.  An  ND2200  4096  channel  analyzer  was 
used  to  obtain  2048  channel  spectra  and  record  the  results 
on  magnetic  tape.  The  tape  was  then  taken  to  the  central 
computer  facility  and  stored  there  for  use  in  computer 
analysis. 

Assay  Equipment 

All  of  the  radioactive  samples  prepared  were  assayed 
in  an  NMC  2rr  geometry  gas  flow  proportional  counter.  The 
associated  electronic  equipment  included  a Hewlett-Packard 
5554A  preamplifier,  an  NIM  Standard  high  voltage  power 


supply,  a Tennelec  TC203BLR  linear  amplifier,  a Tennelec 
TC444  single  channel  analyzer,  and  a Canberra  1772 
counter/ timer. 

A 8 say  Procedure 

The  sample  to  be  assayed  was  inserted  into  the  flow 
counter,  the  counter  sealed,  and  the  air  flushed  out  by  the 
flowing  gas.  Each  sample  was  then  counted  at  least  five 
times  for  one  minute  at  each  of  two  voltages,  1000V  to 
count  just  alpha  particles  and  1600  or  1800V  to  count  both 
alpha  and  beta  particles. 

Preparation  of  Sources 

Kirby  Method.  Most  of  the  sources  prepared  for  this 
study  were  prepared  by  a method  developed  by  H.  W.  Kirby  of 
the  Mound  Facility  in  Miamisburg,  Ohio  (Ref  4 » 134-135).  The 
backing  plate  is  rinsed  for  a few  minutes  in  alcoholic  KOH, 
then  rinsed  with  water  and  dried.  The  area  for  the  sample 
is  masked  with  a circular  rod  while  the  plate  is  sprayed 
with  Krylon.  The  plate  is  then  warmed  and  dried  under  a 
heat  lamp  for  a few  minutes.  The  uncoated  area  is  covered 
with  IN  HNO^  and  allowed  to  stand  for  two  minutes.  The 
plate  is  then  removed,  rinsed  thoroughly  with  distilled 
water,  and  dried  again  under  the  heat  lamp. 

The  radioactive  solution  to  be  deposited  is  dropped  on 
the  uncoated  area  from  a microsyringe  and  allowed  to  evap- 
orate to  dryness.  This  is  then  covered  with  two  drops  of 
0.1  N HNO^  and  heated  for  two  minutes  without  evaporating. 


One  drop  of  3N  NH^OH  Is  then  added  and  the  solution  is 
evaporated  to  dryness.  Two  or  three  sore  drops  of  NH^OH 
are  added  and  again  evaporated  to  dryness.  The  NH^OH  is 
then  driven  off  by  further  heating  on  a hot  plate  until  the 
visible  fumes  are  observed. 

The  major  change  from  Kirby's  procedure  was  to  use  a 
heat  lamp  for  the  initial  drying  operations.  This  was  done 
because  the  available  hot  plates  were  not  able  to  maintain 
a steady  temperature  in  the  90-l00*C  region.  However,  eight 
of  the  first  nine  samples  were  prepared  using  hot  plates. 

Armani  Electroplating  Method.  One  sample  was  prepared 
using  the  "high  voltage"  electroplating  procedure  suggested 
by  Roland  Armani  of  Argonne  National  Laboratory  (Ref  1). 

One  hundred  microliters  of  the  New  England  Nuclear  Radium 
DEF  solution  lot  #72-187,  about  0.1  microcurie,  were  added 
to  10  milliliters  of  2-propanol  in  a 14  millimeter  inner 
diameter  glass  tube,  which  was  clamped  with  rubber  washers 
at  either  end  to  a plate  cleaned  as  for  the  Kirby  method. 

Platinum  gauze  was  not  available,  so  a roughly 
octagonal  plate  of  platinum  about  six  millimeters  in  diameter 
was  used  for  the  anode.  The  anode  was  immersed  just  below 
the  surface  of  the  plating  solution,  which  was  then  plated 
for  two  hours  at  600  volts  and  1.1  microamps  current. 

Only  about  30%  of  the  Radium  DEF  was  plated  out  and  the 
resolution  was  no  better  than  that  of  the  evaporated  samples. 
It  was  decided  not  to  attempt  to  perfect  this  method  and  the 
remaining  samples  were  all  made  by  the  Kirby  method. 
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IV.  Experimental  Results 


Introduction 

This  chapter  covers  the  Ra  DEF  samples  prepared  as 
reference  sources,  the  results  of  the  assay  to  determine 
their  activity,  the  spectra  obtained  from  those  sources, 
and  the  relative  sharpness  of  the  peaks  obtained  from 
those  sources. 

Sources  Prepared 

Twenty-one  sources  were  prepared,  assayed,  and  counted 
with  the  surface  barrier  detector.  Sources  1 and  2 were 
prepared  to  test  the  Kirby  method.  Source  3 was  prepared 
by  electroplating.  It  gave  the  broadest  peak  of  any  of  the 
sources  and  only  20%  of  the  Ra  DEF  in  the  plating  solution 
was  plated  onto  the  planchet,  so  the  remaining  sources  were 
prepared  by  the  Kirby  method.  Two  sets  of  Ra  DEF  samples 
were  then  prepared  with  varying  amounts  of  Pb(N0p2  added  to 
increase  self-absorption. 

Samples  4-9  were  the  first  set  of  samples  of  varying 
thickness.  The  activities  of  these  samples  agreed  within 
±10%.  However,  as  can  be  seen  in  Table  II,  the  spread  in 
the  resulting  spectra  did  not  correspond  to  the  quantity  of 
added  lead  nitrate.  Samples  7-9  were  dried  too  long  on  the 
hot  plate  and  this  apparently  caused  considerable  oxidation 
from  the  nitrate  to  the  oxide.  The  spectra  from  samples  7-9 
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Table  II 


Comparison  of  Full  Width  at  Half  Maximum  In 
2048  Channel  Biased  Spectra  From  100  mm^ 
Surface  Barrier  Detector  with  Quantity 
Added  PbCNO^^  in  the  Sample 


Sample 

# 

Added  Pb(N03)2 
x . 525|Ag  + 37. 

KeV  FWHM 
Unoxtdized 

KeV  FWHM 
Oxidized 

a Act  (Ox) 

a Ac t( Unox) 

1 

0 

39.1 

2 

0 

35.2 

3 

0 

51.6 

4 

0 

39.6 

5 

1 

52.0 

6 

2 

44.6 

7 

3 

24.3 

8 

4 

9 

5 

10 

0 

39.3 

11 

0 

30.2 

12 

1 

29.0 

13 

2 

37.1 

14 

3 

32.8 

15 

4 

35.4 

16 

5 

32.9 

17 

5 

35.6 

18 

0 

39.4 

19 

0 

31.9 

20 

5 

34.2 

21 

5 

37.1 

! 


20.4 

28.1 

34.7 

.90 

29.2 

.60 

35.3 

.89 

34.3 

.87 

34.9 

.92 

35.9 

.98 

32.8 

.87 

had  the  sharpest  peaks  of  all,  but  they  had  double  peaks 
due  to  ringing  in  the  line  from  the  amplifier  to  the 
multi-channel  analyzer  (This  was  corrected  by  Inserting  a 
terminator.)  and  were  not  used  for  computer  analysis. 

Samples  10-21  were  the  second  set  prepared  with  varying 
quantities  of  added  lead  nitrate.  These  were  prepared 
simultaneously  and  dried  under  a heat  lamp  instead  of  on  a 
hot  plate.  The  lead  nitrate  was  added  In  solutions  of 
varying  concentration  so  that  the  same  volume  of  water  was 
evaporated  from  each  of  the  samples.  Samples  1-17  were 
prepared  on  stainless  steel  planchets,  while  18-21  were 
prepared  on  small  squares  of  platinum.  Sources  10-16  were 
later  heated  over  a Fisher  burner  to  oxidize  the  deposit 
from  the  nitrate  to  the  oxide.  These  sources  were  used  for 
the  computer  analyses. 


Assay  Results 

All  activities  are  calculated  based  on  the  NIC  Radium 
DEF  sample  #A110,  with  one  microcurie  equal  to  2.22  x 10^ 
counts  per  minute.  Sample  AllO  was  measured  at  514  disinte 
grations  per  second  on  28  September  1961.  The  activity  on 
24  August  1978  is  computed  from  the  exponential  decay  law. 


A(t) 


c 1/2 


t-t 


A(tQ)  x e 


-(.693)(t-t0) 


cl/2 

20.4  x 365  1/4  days 
330  + 16  x 365  1/4  days 


2; 


A(tQ)  **  514  dps  + 5% 

A( t)  = 289  dps  + 5% 

The  remaining  activities  are  calculated  from  the 


formula, 


Activity  (Sample)  = Counts/mia^Sa^ple),  x Activity  (Standard^ 


As  seen  in  Tables  III  and  IV,  which  give  the  averages 
of  five  one-minute  counts  plus  or  minus  the  standard 
deviation  of  those  counts,  the  unoxidized  samples  agreed 
fairly  well  in  activity  from  sample  to  sample.  As  seen  in 
Table  IV,  the  standard  deviations  of  the  averages  of 
activities  for  samples  on  the  same  backing  (stainless  steel 
or  platinum)  were  only  1-2%.  However,  the  combined  a + 0 
activity  divided  by  three  for  the  three  nuclides  in  the 
decay  chain  was  always  several  percent  higher  than  the  a 
only  activity. 

The  results  of  the  assay  of  samples  10-16  after  oxidizing 
are  listed  in  Table  V.  The  combined  a + 0 activity  remained 
consistent  from  sample  to  sample.  However,  the  fluctuation 
In  a only  counts  increased  to  a standard  deviation  of  15% 
as  varying  amounts  of  polonium  were  driven  off. 


Spectral  Results 

The  widths  of  the  spectra  from  all  the  samples  are 
compiled  in  Table  II.  There  was  no  correlation  between  the 
amount  of  added  lead  nitrate  in  a sample  and  the  full  width 


at  half  maximum  of  the  peak  from  that  sample.  Most  of  these 
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Table  III 


Results  of  Assay  of  Samples  1-9  in  2-rr  Geometry 
Gas  Flow  Proportional  Counter  at 
1000V  (a)  and  1600V  (a  + g) 


Sample  a Counts/60  sec  a + 3 cts/60  sec  Activity 


a (a+9 )/3 


1 

78,490 

+ 

170 

295,300 

+ 

2100 

.0721 

.0846 

2 

77,800 

+ 

230 

336,920 

± 

640 

.0715 

.0966 

3 

21,730 

+ 

170 

112,990 

± 

290 

.0200 

.0324 

4 

122,640 

+ 

390 

470,150 

± 

780 

.113 

.135 

5 

122,570 

± 

460 

466,580 

± 

510 

.113 

.134 

6 

120,520 

± 

300 

474,200 

± 

1300 

.111 

.136 

7 

102,030 

± 

250 

455,450 

± 

470 

.0938 

.131 

8 

123,960 

± 

220 

483,180 

± 

480 

.114 

.138 

9 

99,500 

± 

300 

402,830 

± 

680 

.0914 

.115 

AllO 

8,512 

± 

92 

27,291 

± 

165 

.00782 

.00782 

+.00039 

NOTE i Samples  1-9  were  prepared  on  stainless  stell 

planchets.  Sample  3 was  prepared  by  electroplating. 
Sample  AllO  was  used  as  a calibration  standard. 


Table  IV 

Results  of  Assay  of  Samples  10-21  in  2v  Geometry 
Gas  Flow  Proportional  Counter  at 
1000V  (a)  and  1800V  (a  +3) 


Sample  a Counts  a + 9 Counts  a Act  (a+9)/3  Act 

# per  60  seconds  per  60  seconds  pC^ 


10 

109,780 

+ 

180 

363,120 

+ 

560 

.101 

.104 

11 

107,470 

+ 

230 

359,490 

+ 

800 

.0988 

.103 

12 

109,310 

+ 

200 

373,910 

+ 

660 

.100 

.107 

13 

111,740 

+ 

170 

377,640 

+ 

560 

.103 

.108 

14 

108,740 

+ 

310 

371,480 

+ 

520 

.0999 

.106 

15 

110,770 

± 

320 

369,880 

+ 

670 

.102 

.106 

16 

108,450 

± 

280 

352,250 

+ 

850 

.0997 

.101 

17 

110,860 

± 

140 

363,600 

+ 

6300 

.102 

.104 

18(pl) 

116,210 

+ 

150 

393,780 

+ 

380 

.107 

.113 

19(p2) 

113,410 

+ 

450 

390,340 

+ 

490 

.104 

.112 

20(p3) 

115,090 

+ 

320 

403,070 

+ 

470 

.106 

.116 

21(p4) 

113,730 

+ 

240 

398,810 

+ 

650 

.105 

.114 

Average 


--  - — — — o — 

10-17 

.1008 

.1049 

+.0014 

+.0023 

Average 

18-21 

.1055 

.1138 

+.0013 

+.0017 

NOTEi  Samples  10-17  were  prepared  on  1.25-inch  diameter 

stainless  steel  planchets,  while  18-21  were  prepared 
on  14  mm  square  platinum  planchets,  which  were  then 
placed  on  stainless  steel  planchets  for  counting. 


30 


Table  V 


Results  of  Assay  of  Samples  10-16  In  2it  Geometry 
Gas  Flow  Proportional  Counter  at 
1000V  (a)  and  1800V  (a  + 0) 


Sample  a Counts 

# per  60  seconds 

a + 0 Counts 
per  60  seconds 

Activity, 

a 

fiC  t 

Activity, 

(a+0)/3 

10 

99,080 

+ 

460 

381,680 

+ 

120 

.0911 

.109 

11 

64,860 

+ 

220 

342,390 

+ 

180 

.0596 

.0981 

12 

97,070 

+ 

240 

367,080 

+ 

540 

.0892 

.105 

13 

97,300 

± 

310 

381,660 

610 

.0894 

.109 

14 

99,860 

+ 

200 

361,360 

510 

.0918 

.104 

15 

108,980 

+ 

310 

371,930 

+ 

790 

.100 

.107 

16 

94,260 

+ 

590 

368,150 

+ 480 

.0866 

.106 

Average 

.087  + 
.013  “ 

.1054  + 
.0037  ~ 

NOTE i These  samples  were  gently  heated  over  a Fisher 
burner  to  convert  the  deposit  from  the  nitrate 
to  the  oxide. 
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peaks  were  10-15  KeV  wider  than  the  21  KeV  wide  peaks  from 
the  Ani  1 and  Pu  #2  samples,  which  are  plotted  in 
Figures  6 and  7 . However,  there  was  a good  selection  of 
widths  in  the  region  of  #1.  Its  full  width  at  half 

maximum  was  about  32  KeV,  but  could  not  be  determined 
precisely  because  the  notch  between  the  twc  main  peaks  did 
not  come  down  to  one-half  the  maximum  (See  Figure  7), 

Because  of  this  good  match,  Pu'1  #1  was  selected  as  the 
"unknown"  to  be  analyzed  by  the  computer  program  using  as 
reference  peaks  the  spectra  from  the  oxidized  sources  10-16, 
which  are  plotted  in  Figures  8-15. 
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AM  241  ORTEC  1/21/69  0.6  KEV/CHANNEL  1C. 000  SEC 


FIG-  7.  PU  240  * 2 flRMRNl  77  0.6  KEV/CHRNNEL  10.000  SEC. 


FIG  8.  PU  240  * 1 ARMAN l 77  0.6  KEV/CHANNEL  20.000  SEC 
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FIG-  14.  Rfl  OEF  * IS  (OXIDIZED)  0.6  KEV/CHANNEL  40.000  5EC 


FIG-  15.  RR  DEF  • 16  (OXIDIZED)  0.6  KEV/CHRNNEL  40,000  SEC 


V . Program  Modifications 


Introduction 

The  starting  point  for  these  programs  was  Program 
ALPHAFIT,  which  was  produced  by  Lieutenant  Richard  S.  Hartley 
with  the  help  of  Phil  Poirier  from  Program  AUTOFIT  written  by 
J.  R.  Comfort  (Ref  5).  Two  specific  modifications  were 
originally  planned,  to  convert  the  program  to  read  the 
spectral  data  directly  from  magnetic  tape  or  disc  storage, 
and  to  use  a library  of  curves  to  choose  the  curve  which 
will  give  the  best  analysis  and  use  it  to  analyze  the  unknown 
spectrum.  In  addition,  subroutines  were  written  to  create 
Gaussian  spectra  modified  by  a Poisson  distribution  for 
counts  in  each  channel  and  these  spectra  were  used  in  place 
of  detector  spectra  to  study  the  program's  performance. 

Conversion  to  Tape 

It  was  desired  to  convert  the  program  to  read  the 
spectral  data  directly  from  the  magnetic  tape  on  which  they 
were  written  by  the  multi-channel  analyzer.  Two  programs, 
SCRIPT  and  SPECTRA  were  written  by  H.  Careway  (Ref  3)  to  use 
the  FORTRAN  BUFFER  IN  and  DECODE  statements  to  read  these 
data  into  an  array.  SCRIPT  then  prints  out  the  numbers, 
while  SPECTRA  produces  a CALCOMP  plot.  The  READ  section  of 
these  programs  was  extracted  and  modified  into  subroutine 
RDSPEC.  RDSPEC  was  then  called  directly  from  the  main 
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program  to  read  the  reference  spectra  and  from  subroutine 
SPCTRM  to  read  the  unknown  spectrum. 

Reading  the  data  in  this  manner  La  not  without  diffi- 
culty. First,  the  spectra  to  be  read  must  be  placed  in 
order  on  file  TAPE  1 by  using  the  SK.IPF  and  COPYCF  control 
card  statements.  Then,  the  statement  FILE,  TAPEl,  FO=i>Q, 
BT'C,  RT=S,  BFS=1000.  is  used  to  tell  the  computer  this  is 
a stranger  tape  (RT=S)  and  the  buffer  size  (BFS=1000).  This 
must  be  followed  by  a LOAD-SET  statement  such  as  LDSET, 

FILES  =TAPE 1 , PRESET  =ZERO . 

Second,  perhaps  because  of  a stray  bit  or  an  unclean 
tape,  a datura  would  be  rejected  by  the  computer  dumping  the 
program  immediately.  RDSPEC  was  modified  with  the  statement, 
ERRSET  (FOUNTS, 100) , so  that  the  program  would  not  dump  for 
just  one  such  datum.  The  numbers  read  in  this  way  were 
always  satisfactory,  but  occasionally  the  bottom  section  of 
a spectrum  would  be  rejected  by  the  computer. 

Reorganization 

Program  ALPHAFIT  analyzed  spectra  of  2048  channels  for 
up  to  20  separate  peaks.  This  is  necessary  for  analyzing 
a spectrum  with  several  isotopes  over  a wide  range  of 
energies  but  requires  an  enormous  (over  250k)  quantity  of 
core  storage.  This  requirement  impedes  studying  the 
behavior  of  the  program,  making  changes,  and  troubleshooting, 
so  the  program  was  reduced  to  analyzing  512  channels  out  of 
the  2048  containing  up  to  five  peaks.  This  reduced  the 
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core  requirements  to  120-124K,  which  allowed  several  runs 
each  evening.  However,  this  reorganization  imix>sed  certain 
limits.  With  the  biased  amplifier  spreading  the  spectrum 
out  to  0.48  or  O.bO  ReV  per  channel,  a range  of  only 

241 

250-300  KeV  could  be  studied.  This  was  adequate  tor  am 

240 

but  insufficient  for  the  Pu  spectra  taken,  so  only  the 

240 

two  largest  peaks  in  Pu  were  analyzed.  More  than  five 
peaks  would  have  to  be  analyzed  in  most  unknown  spectra. 

The  second  change  was  to  have  the  program  read  and 
store  five  reference  spectra  to  analyze  each  unknown 
spectrum.  This  was  done  simply  by  copying  each  spectrum 
into  the  array  REFSTO  (520,5)  before  using  the  reference  to 
analyze  the  unknown. 

Artificial  Spectra 

In  order  to  study  the  performance  of  the  program's 
analyzing  routines,  it  was  desired  to  generate  spectra  of 
arbitrary  characteristics  to  simulate  alpha  spectra.  Sub- 
routine CONJUR  was  written  to  use  several  routines  written 
by  Bevington  (Ref  2)  to  generate  Gaussian  peaks  of  speci- 
fied height,  standard  deviation  and  position,  add  them  to 
create  a spectrum,  then  randomly  modify  the  counts  in  each 
channel  according  to  an  approximate  Poisson  distribution. 
CONJUR  and  the  associated  subroutines  were  then  substituted 
for  subroutine  RDSPEC.  The  resulting  program  was  called 
ALFAIC. 


VI. 


Program  Computation  Results 


Introduction 

The  program's  analyzing  performance  was  extensively 
tested  under  idealized  conditions  by  using  generated 
Gaussian  type  spectra.  Samples  of  americium  and  plutonium 
were  analyzed  by  the  512  channel  version  using  the  best  of 
the  samples  prepared  (after  oxidizing).  One  of  the  pluto- 
nium samples  was  also  analyzed  using  the  2048  channel 
version  of  the  program.  These  analyses  showed  that  a 
difference  of  only  one  channel  between  the  widths  of  the 
reference  and  unknown  peaks  caused  significant  reduction  in 
the  accuracy  of  the  analysis. 

Conjured  Spectra 

The  program  was  used  to  analyze  spectra  composed  of 
two  Gaussian  peaks  with  a standard  deviation  (o)  of  20 
channels.  The  separation  between  these  peaks  was  varied 
from  one  quarter  to  two  and  a half  times  the  standard 
deviation  (5-50  channels)  to  investigate  the  anomaly 
reported  by  Hartley  (Ref  5i22).  The  width  of  the  reference 
peak  was  varied  from  15  to  25  channels  to  observe  this 
effect. 

Separation  Effect.  The  program  was  tested  with  peaks 
of  5 » 10  and  1 < 10  relative  heights  (Tables  VI  and  VII).  None 
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Table  VI 


Errors  Versus  Separation  Distance  for  Gaussian 
Peaks  1/2  = 5il0,  a = 20,  Peak  1 in 
Channel  250,  True  Area  ~ 572000 


Peak  2 


Area( true )-Area( comp) 

RMS 

Estimated 

Area (true) 

7. 

Error  % 

Peak 

Error 

1,  Peak  2 ) 

.064 

.035 

805 

865 

.065 

.038 

525 

543 

.065 

.037 

355 

390 

.062 

.051 

313 

357 

.0604 

.046 

254 

301 

.0581 

.035 

196 

236 

.0538 

.033 

164 

209 

.0480 

.035 

152. 

5 202.4 

Ratio 
Peak  2 


of  these  peaks  were  sufficiently  separated  to  show  a drop 
in  counts  between  the  two  peaks  (Figures  16-17).  (Compare 
with  the  americium  and  plutonium  spectra  in  Figures  6 
and  8 . ) 

The  absolute  difference  between  the  computed  and  true 
areas  rises  gradually  as  the  two  peaks  approach.  When  the 
peaks  are  separated  by  less  than  one  standard  deviation,  the 
Root  Mean  Square  difference  between  the  computed  approxi- 
mation and  the  analyzed  spectrum  declines  as  the  smaller 
peak  vanishes  into  the  larger,  while  the  computed  error  in 
each  peak  is  generally  larger  and  the  computed  ratio  of  the 
two  peaks  fluctuates  wildly. 

These  results  show  the  program  is  unable  to  resolve 
two  peaks  accurately  when  they  are  so  close  as  to  appear  to 
be  one  (See  Figure  18).  Nevertheless,  both  the  difference 
between  the  absolute  areas  of  the  unknown  and  resolved 
spectra,  the  absolute  residual,  and  the  Root  Mean  Square  of 
the  difference  at  each  point,  the  average  residual,  remain 
less  than  0.1%.  The  two  peaks  have  merged  and  are  indistin- 
guishable. The  only  indication  that  the  spectrum  is 
composed  of  more  than  one  peak  is  the  1-10%  difference 
between  the  left  and  right  half  widths  at  half  maximum. 

Reference  Width  Effect.  The  analysis  of  peaks  of  fixed 
width  with  references  of  various  widths  show  (See  Tables  VIII 
and  IX  and  Figures  19-20)  that  matching  the  spread  of  the 
reference  to  that  of  the  unknown  can  reduce  the  residuals 


Fig.  16.  Randomized  Gaussian  Peaks.  Separation  ■ 50  Channels 
Heights  are  1000,  10000. 


Fig.  18.  Randomized  Gaussian  Peaks.  Separation  » 20  Channels 
Heights  are  5000,  10000. 


Table  VIII 


Analysis  of  (a  = 20)  Gaussian  Peaks  with  Reference 
Peaks  of  Varying  a.  Ratio  of  Peak 
Heights  = 10 »1  Separation  = 50(2.5  o) 


Ref  erence 
a 

True  Area 
-Resolved  Area 

RMS 

Error 

Resolved  Ratio 
of  Two  Peaks 

15 

2092 

10299 

7.1 

16 

1175 

7857 

7.6 

18 

134 

3570 

8.6 

19 

296 

1374 

7.9 

19.5 

4.2 

703 

8.8 

20.0 

237 

171 

10.05 

20.5 

353 

761 

11.1 

21.0 

436 

1492 

12.3 

22 

564 

2954 

13.5 

24 

731 

5363 

18.3 

25 

703 

11297 

-66 

Table  IX 


Analysis  of  (a  = 20)  Gaussian  Peaks  with  Reference 
Peaks  of  Varying  a.  Ratio  of  Peak 
Heights  = 10 i 5 Separation  = 50(2.5  a) 


— — 

Reference 

True  Area 

RMS 

Resolved  Ratio 

0 

-Resolved  Area 

Error 

of  Two  Peaks 

15 

4392 

13127 

1.72 

1 

16 

2527 

10130 

1.77 

18 

623 

4724 

1.88 

19 

6.4 

1926 

1.939 

19.5 

168.6 

964.8 

1.9699 

20.0 

272.8 

198.4 

2.0008 

20.5 

394.6 

927.8 

2.046 

21.0 

500.2 

1728.3 

2.11 

22 

678 

4154 

2.23 

24 

952 

7856 

2.56 

25 

1072 

9580 

2.76 

54 


Standard  Deviation  of  Reference  Peak. 

Fig,  19.  Analysis  of  (a  = 20)  Gaussian  Peaks  with  Reference 
Peaks  of  Different  Widths.  Ratio  of  Peak  Heights  = 10 i 1 . 


Standard  Deviation  of  Reference  Peak 

Fig.  20.  Analysis  with  Reference  Peaks  of  Different  Widths 
Standard  Deviation  of  Peaks  Analyzed  = 20.  Ratio  of  Peak 
Heights  = 10 i 5. 


by  two  or  three  orders  of  magnitude  under  ideal  conditions. 

A difference  of  only  0.5  in  the  standard  deviation  produces 
a clear  increase  in  the  RMS  error. 

Real  Spectra 

512  Channel  Analyses.  The  first  group  of  analyses  were 
conducted  on  a 512  channel  segment  out  of  the  full  2048 
channel  spectrum.  Five  spectra  were  used  to  analyze  two 
samples,  Am*^*  and  Pu^^  #1.  The  spectra  from  the  oxidized 
Ra  DEF  samples  10,  11,  12,  13,  and  16  were  selected  because 
they  had  the  sharpest  peaks.  Although  these  peaks  could  not 
be  placed  in  sequence  showing  a continuous  regular  change 
in  the  curve  shape,  they  could  be  placed  in  order  of  their 
full  width  at  half  maximum.  As  shown  in  Table  X and 
Figures  21  and  22,  the  RMS  error  corresponded  directly  to 
the  fit  between  the  FWHM’s  of  the  unknown  and  the  reference 
used  for  analysis.  It  first  dropped  and  then  rose  as  the 
FWHM  increased  for  the  plutonium  sample  and  rose  steadily 
for  the  americium  sample. 

Other  measurements  of  the  accuracy  of  the  fit  were  not 

as  well  correlated  to  the  FWHM.  The  difference  in  areas 

when  analyzing  Pu*-  #1  showed  the  same  trend  as  the  RMS 

error  except  for  sample  #10,  but  there  was  no  agreement 

when  analyzing  americium.  The  computed  ratio  of  the  two 
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main  peaks  in  Pu  was  too  low  except  for  sample  16. 

2048  Channel  Analyses.  The  analyses  on  the  full  2048 

240 

channel  spectrum  of  Pu  showed  good  correlation  between 
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Table  X 

Analysis  of  Plutonium  and  Americium  with  Polonium 
Spectra  of  Varying  Self -Absorption.  All 
Spectra  from  100  mra^  ORTEC  Surface  Barrier 
Detector  at  30-35  V Net  Bias 


Unknown  = Pu^^  vl 

FWHM  = 

about  32 

Reference  Spectrum 
from  Sample  # 

FWHM 

(KeV) 

A Area 
Error 
(counts) 

RMS 

Error 

(counts) 

Computed 

Ratio 

11 

29.2 

9318 

11244 

2.79 

16 

32.8 

7833 

10618 

3.06 

13 

34.3 

7206 

6943 

2.42 

10 

3A.7 

23038 

5942 

2.67 

12 

35.3 

7773 

9616 

2.71 

?4 1 

Unknown  = Am 

FWHM  = 21. 1 

Reference  Spectrum 
from  Sample  if 

FWHM 

(KeV) 

A Area 
Error 
(counts) 

RMS 

Error 

(counts) 

11 

29.2 

29140 

18325 

16 

32.8 

24330 

20954 

13 

34.3 

15630 

21835 

10 

34.7 

9642 

21709 

12 

35.3 

19826 

23368 

NOTE i Only  512  channels  of  these  spectra  containing  the 
two  or  three  chief  peaks  were  analyzed.  The  ratio 
for  plutonium  is  for  the  5.168  and  5,123  meV  peaks. 
The  true  ratio  is  3.167. 
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RMS 


FWHM(KeV)  of  Reference  Spectrum 

Analysis  of  Pu^O  #1  with  Various  Ra  DEF  Samples 
•a  From  ORTEC  100  mm^  SBD  at  30-35  V Net  Bias. 


RMS 


FWHM(K.eV)  of  Reference  Spectrum 

Fig.  22.  Analysis  of  Am241  with  Various  Ra  DEF  Sample 
All  Spectra  From  ORTEC  100  mm2  sqD  at  30-35  V Net  Bias 


all  three  measures  of  the  fit  and  the  match  of  the  FWHM's 
except  for  reference  sample  #15  (See  Table  XI  and 
Figure  23).  This  had  the  largest  FWHM  but  gave  the  second 
best  fit.  The  resolved  composite  areas  were  50%  to  150% 
too  large  and  the  composite  now  gave  too  large  a percentage 
of  the  counts  to  the  largest  peak.  These  results  confirm 
that  matching  the  reference  curve  shape  to  the  spectrum  being 
analyzed  gives  a more  accurate  analysis. 

The  results  of  these  analyses  show  the  potential  of  the 
library  of  curves  approach  to  produce  a more  accurate 
analysis  of  alpha  spectra.  However,  the  program  remains 
unstable  and  capable  of  producing  bizarre  results  even  under 
good  conditions.  The  peak  locating  routines  will  sometimes 
shift  the  locations  of  two  clearly  separated  peaks  to  only  a 
few  channels  apart  and  at  other  times  will  shift  one  or 
more  peaks  completely  out  of  the  spectrum.  The  routines 
which  compute  the  relative  multipliers  for  each  peak  (See 
Eq  (4))  will  at  times  return  a negative  value  for  one  of 
the  peaks.  These  problems  need  to  be  solved  in  future  work 
with  this  program  before  it  can  be  considered  completely 
reliable.  Even  then,  the  program's  performance  will  be 
limited  by  the  reference  spectra  available. 


Table  XI 


Results  of  2048  Channel  Analyses  of  Pu24^  #1 


Pu240  FWHM  * 103.58  Channels 


Reference 

Sample 

FWHM 

Channel  s/KeV 

A Area 
Counts  /7» 

RMS 

Count  s/7, 

Peak 

Ratio 

llOx 

50.68/29.2 

1624448/142 

75965/6.6 

16.4 

160x 

57.07/32.8 

1033675/90 

53577/4.7 

7.3 

130x 

59.59/34.3 

551953/48 

37264/3.25 

4.42 

lOOx 

58.05/34.7 

902232/79 

46422/4.05 

5.8 

120x 

60.59/35.3 

1203889/105 

63790/5.6 

17.4 

140x 

60.63/35.3 

1375767/120 

63419/5.5 

12.4 

150x 

62.47/36.4 

831737/73 

43413/3.8 

4.98 

NOTE i The  peak  ratio  is  for  the  5.168  and  5.123  meV  peaks. 
The  true  ratio  is  3.167. 


Reference  FWHM  (KeV) 

240 

Fig.  23.  Results  of  2048  Channel  Analyses  of  Pu 
with  Various  Ra  DEF  Samples.  All  Spectra  From  ORTEC 
100  mm2  SBD  at  30-35  V Net  Bias. 
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VII . Conclusions  and  Recommendations 


Conclusions 

The  analyses  with  generated  spectra  clearly  show  that, 
at  least  under  Ideal  conditions,  matching  the  FWHM  of  the 
reference  peak  to  the  FWHM  of  the  spectrum  to  be  analyzed 
gives  a better  analysis.  Differences  of  only  one  channel 
In  these  FWHM's  were  shown  to  Increase  the  RMS  error  by  a 
factor  of  3 to  5.  This  strongly  supports  the  need  for  a 
library  of  reference  curves  to  use  in  analyzing  various 
alpha  spectra. 

The  analyses  with  real  spectra  also  support  these  con- 
clusions. The  reduced  range  (512  channel)  analyses  of 
americium  and  plutonium  both  show  a reduction  in  the  RMS 
error  as  the  FWHM  of  the  reference  peak  is  made  closer  to 
the  FWHM  of  the  spectrum  analyzed.  The  full  range  (2048 
channel)  analyses  of  plutonium  confirm  this  reduction  of 
the  RMS  error  and,  in  addition,  show  a corresponding 
reduction  in  the  absolute  residual  area  and  more  accurate 
resolution  of  the  relative  areas  in  the  two  chief  peaks. 

These  analyses  show  that  the  use  of  a library  of 
reference  spectra  can  give  a several  percent  reduction  in 
error  provided  that  reference  spectra  with  the  correct 
characteristics  are  available.  The  chief  need  at  this  time 
is  for  a consistent  set  of  smoothly  varying  spectra, 


including  a range  of  FWHM's  as  low  as  20  KeV  or  better. 
Recommendations 

Samples.  The  main  need  is  for  better  reference  samples. 

210 

Po  is  clearly  the  isotope  to  use.  It  may  be  possible  to 
prepare  adequate  Ra  DEF  samples,  but  these  will  have  to  be 
the  oxide  or  perhaps  bare  metal.  The  nitrate  samples  have 
too  much  self -absorption,  even  with  no  added  lead.  An 
electroplating  technique  will  probably  be  required  to 
produce  samples  with  low  enough  self -absorption. 

There  are  two  further  studies  that  could  be  made  with 

the  samples  prepared  for  this  thesis.  New  spectra  should 

be  obtained  frort  samples  7-9.  If  the  peaks  from  these  have 

FWHM's  as  narrow  as  indicated  in  Table  II,  they  should  be 

used  to  analyze  samples  with  correspondingly  low  self- 

absorption  such  as  the  americium sample  used  in  this 

study.  Additional  spectra  from  the  later  set  (10-21)  should 

also  be  taken  in  another  year  or  two  to  check  for  changes  in 

the  curve  shapes  due  to  a more  even  distribution  in  the 
210 

polonium  as  a result  of  radioactive  decay. 

Program.  There  are  also  improvements  that  can  be  made 
in  the  program.  A big  problem  is  how  to  store  and  access 
a large  library  of  2048  channel  reference  spectra.  These 
will  probably  have  to  be  stored  on  magnetic  tape,  so  some 
means  of  conveniently  reading  a particular  record  is 
needed. 


The  operation  of  the  main  program  could  be  considerably 
streamlined  by  eliminating  some  of  the  unused  options.  The 
following  sequence  should  be  adequate.  First,  controlling 
data  are  read.  Second,  the  spectrum  to  be  analyzed  and  a 
background  spectrum  are  read.  The  FWHM  is  computed  for  the 
largest  peak  in  the  unknown  spectrum  and  a reference 
spectrum  is  selected  and  read.  The  program  then  passes  these 
three  arrays  to  the  analysis  routines  and  conducts  the 
analysis.  The  automatic  and  interpolated  background  options 
and  the  computer-derived  peak  shape  options  could  then  be 
eliminated. 

It  would  also  be  useful  to  modify  the  analyzing  routines 
so  that  the  spectrum  for  each  peak  extends  all  the  way  to 
zero,  perhaps  using  a computed  average  of  the  lower  channels 
to  fill  the  gap  from  the  bottom  channel  down  to  zero. 

Finally,  the  analysis  routines  should  be  carefully  studied 
to  verify  their  correct  operation  and  make  useful  modifi- 
cations such  as  rejecting  negative  peaks  and  limiting  the 
shifting  of  the  specified  peaks  to  a fixed  range  about  the 
specified  location.  These  routines  are  the  subroutines 
PREPAR  and  DAVIDN  and  all  the  subroutines  called  by  these 
two. 
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Appendix  A 


Utility  Programs 


Introduction 

Several  programs  besides  ALPHFT  (alpha  fit)  were  used 
during  the  work  on  this  thesis.  It  was  decided  to  read  the 
multi-channel  analyzer  generated  data  from  magnetic  tape 
(or  tape-derived  permanent  disc)  storage.  Captain  Harold 
A.  Careway  (Ref  3)  had  written  two  programs,  SCRIPT  and 
SPECTRA,  to  perform  this  operation  using  the  FORTRAN  BUFFER 
IN  and  DECODE  statements. 

Program  SCRIPT 

Program  SCRIPT  was  designed  to  read  in  a series  of  L 
consecutive  spectra,  then  print  out  the  first  IC  channels  of 
each  spectrum.  In  the  form  developed  by  Careway,  both 
SCRIPT  and  SPECTRA  were  quite  sensitive  to  some  phenomenon 
which  caused  the  code  for  a number  to  be  interpreted  as  the 
code  for  (,  @,  ",  or  some  other  character.  This  would 
cause  the  computer  to  dump  the  program  without  reading, 
writing,  or  plotting  any  more  spectra.  At  the  suggestion  of 
Captain  Stinson,  these  programs  were  modified  using  the 
CALL  ERRSET  statement  to  allow  the  computer  to  accept  a few 
such  events  before  dumping  the  program.  This  allowed  sub- 
sequent spectra  to  be  read  and  written  or  plotted  on  the 
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same  computer  run.  By  printing  out  only  the  first  few  lines 
from  a spectrum,  program  SCRIPT  could  easily  be  used  to 
investigate  which  spectra  contain  bad  data.  However,  when 


the  full  spectrum  was  printed  out  after  such  an  occurrence, 
no  bad  data  were  observed. 

Program  SPECTRE 

Program  SPECTRA,  the  plotting  routine,  was  extensively 
modified  so  that  it  could  be  used  to  plot  the  graphs  for 
this  thesis.  It  was  then  renamed  SPECTRE.  In  the  original 
form,  it  generated  only  a large  graph  of  variable  dimensions 
similar  to  the  type  ?.  plot  in  SPECTRE  but  with  no  title.  The 
program  always  plotted  channels  70  to  IC,  an  input  variable. 
The  spectra  to  be  plotted  had  to  be  placed  in  order  on  the 
file  TAPEl  by  correct  use  of  the  SKIPF,  COPYCF,  and  REWIND 
control  cards.  Provision  was  made  to  edit  erroneous  data 
prior  to  plotting. 

Modifications 

The  program  was  first  modified  using  CALL  ERRSET  as  was 
SCRIPT.  It  was  then  modified  using  additional  CALL  PLOT 
statements  and  a CALL  SYMBOL  statement  to  generate  a smaller 
plot  enclosed  in  a box  with  boundaries  marked  for  cutting 
with  a title  at  the  bottom  of  the  page  (after  cutting).  The 
titles  were  initialized  in  the  array  ITITLE  by  a series  of 
DATA  statements  and  the  correct  title  was  determined  from 
the  tagword  of  the  data  to  be  plotted.  Finally,  the  program 
was  heavily  modified  to  allow  the  user  to  read  the  4XN  array 
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NTAG,  whi.cn  specified  the  spectrum/tagword  to  be  plotted,  a 
choice  of  the  old  and  new  formats,  and  the  first  and  last 
channels  to  be  plotted.  It  is  still  desirable  to  register 
the  first  spectrum  to  be  plotted  by  using  the  SKIPF  control 
card  in  order  to  reduce  Input/Output  time,  but  this  is  not 
necessary.  The  tagwords  to  be  plotted  may  be  read  in  any 
order  but  will  be  plotted  in  ascending  order. 


BE  REGISTERED  BY  USE  O'7 

1=  DUMBER  OP  SPECTRA 
IC=  NUMBER  OF  CHANNELS 

AT.  •»  , « * i *.  .,  i t.  * V » i.  » * ■ tO'f*  t 1 •*•»  + *•** 

RPAO  1 C DC  , L 
IvEAD  1 ">  30  » IE 
DO  1C  n J= 1, L 
N=1 
KK  = 1 

10  BUFFER  IN  (1 ,0)  MX(  1>  ,MW(160)  ) 

*<=Uf  IT (1) 

CALL  EERSET  CROON  i ,10  0) 

PECOrE  (i>  0 , 5 P , )W)H,  M*,X 
IF (KfUNT. GY . 0)  GO  TO  1 ° 0 

SO  FORMAT  (2X,I  i 0 ,9  ( /25r=>.  0)  ,/,f“E».0) 

DO  SO  1 = 1, .?«  & 

Y(N)  = y (i) 

8 0 N = N + 1 

IF(KR.NE.C) GC  TO  10 

PRINT*,"  UNIT  ( 1 ) = ", K^,"  N = ",N 

PRINT  1013,  M 

NN=(?C/10>  +1 

DO  90  ^=1 , NN 

IN= (K-l) ' 10 

T T = If  +1 

PRIN'r  1 0 20,  IN, Y ( T I ) , Y < TI  + 1),Y(II«-2),Y<  II  «G),Y(II+4),  Y(II*5> ,Y  (11  + 
C6J  , Y ( 11+/)  , Y ( 1 1 <•  c ) , Y (11*9) 

DO  CONTINUE 

10"  CONTINUE 
10r’0  FORMAT  (3  3) 

1010  FORMA T < 1 II , TS1 ,"T A 9W3RO  NUM3 ER"» T 55 » [ 5) 

10’C  FORMAT  (1X,I[>,10F10»0) 

STCD 

END 


pROGr  AH  SCRIPT (INPUT, output , TAP-1 ) 

DIMENSION  Y ( 4 1 0 0 ' , X ( ?Zr')  ,*'4  ( lbQT 
C.  this  PROGRAM  prints  out  spectra 
C TMr  FITST  SPECTRA  TO  PRINTED  OUT  UST 
C T J r SRIPP  CONTROL  CARD 

r tjp  pins’  data  card  ts  i"  pormat  ***** 

r T'JP  SECOND  r'  (■  T A CARD  IS  IE  FORMAT  **** 
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Program  SCRIPT  is  simple  to  use.  The  spectra  to  be 


read  are  placed  on  file  TAPEl.  The  FILE  and  LDSET  statements 
tell  the  computer  the  buffer  size  for  the  file  and  load  it  to 
be  read.  Two  data  cards  are  required  in  I 5 format.  The 
first  should  contain  L,  the  number  of  spectra  to  be  read  and 
printed.  The  second  should  give  IC,  the  number  of  channels 
to  be  printed  out  for  each  spectrum  starting  with  channel  0. 

A typical  control  card  sequence  would  be  the  following. 

FTN. 

ATTACH,  TAPEl,  XFILE,  MR=1. 

FILE,  TAPE 1 , F0=SQ,  BT=C,  RT=S,  BFS=1000. 

LDSET,  FILES=TAPEl , PRESET =ZERO. 

LGO. 


1 
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0 0 0 0 0%.  •)  o 


F<(.GC  AM  JPECHECNPUT  , OUTPUT  , PLOT,  TAPED 

V-»  « > 4 *.4  4 » 4 » >»■».*  t » , 

P F U G -d  M SPECT-E  IS  * N ELAP0RATr3  VERSION  OF  PROGRAM  SPECTRA  , 
u <-•  I r m w.\c  provided  "y  np.  GEORSE  john. 

IT  ALLOWS  VH-  =L3TTTNG  OF  TWO  3IFFEdENT  PIO’S  INTERCHANGEABLY! 
ONE  IS  A SMALL  °ORCEkcD  p.OT  IN  AFj.7  TMESIS  FORMAT. 

THE  L 7 HE  T TS  [,  LARGER  V'RSION,  WHICH  “AY  PE  MORE  SUIT  A PL  E 
FOP  Pj;  ELSE  T 30MrARI50NS. 


01  vrf  SION  XI 410  J>  , Y(  1 IPO)  , 
CIlsEFSIO"  IT  I TL  E (b  » 1 1 0 ) 
PlMEfSIGN  NTAS  (-♦,  TO)  , XA  (2) 

(J4*»44»»4;  J 1,44  »■»»»*  444 

ntlf  1TITLE(  1,  H2) /<•  1HFTG 

■>  ( 2tb> 

, YA  ( 2) 

4 4 « 4-  * -* 

. PA 

* w 4 ( 160) 

Mr>,4.«* 

3 E F it  17 

!>♦•*  *****  * * 

( JNOYIDI7ED) 

0.6 

KEV/CHAN 

CNEL  ‘ 0 0 3 C SEC.  / 

DATA  ITITLE(l,s3) /MOHFIG 

• 

RA 

3 E F 

n 

IT 

( J NCXI0I7E0) 

0.6 

KEV/CHAN 

CNEL  ‘ 3L0C  SEC.  / 

OAT/1  ITI  fLEfli  •++)/6IHFjG 

• 

PA 

3" 

« 

13 

( JNPXT0I7E0) 

0.6 

KEV/CHAN 

CNF. L JCCJL  SEC.  / 

DATA  ITITLE( 1,45) /&OHFIG 

• 

RA 

3 E r 

» 

13 

( JNOXIOIZEO) 

0.6 

KEV/CKAN 

CNEL  1CTJC  SEC.  / 

DATA  ITITLEd,  46) /faTHFIG 

• 

RA 

3Ee 

it 

13 

( JNOXIOIZEO) 

0.6 

KEV/CHAN 

CNEL  1UUJ0  SEC.  / 

DATA  ITI T L E ( 1,4?) / G 3 H F T G 

• 

RA 

3r  F 

n 

19 

( JNOXIOIZEO) 

0.6 

KEV/CHAN 

CNEL  1 0 00  0 SEC.  ' . 

CAT;.  ITI  i L E (1 , 3 ) / 6 3 h F ~ g 

• 

RA 

3 E f 

n 

20 

(JNOXIOIZEO) 

0.6 

KEV/CHAN 

CNEL  1 0000  SEC.  / 

DATA  ITITlE(1,h9)/61HFTG 

• 

PA 

3 E F 

it 

20 

( JNOXIDIZEO) 

0.6 

KEV/CHAN 

CNEL  lUOOl  SEC.  ' 

DATA  ITIT Lr(l.bfl) /63HFIG 

• 

R A 

3EF 

it 

21 

( JNOXIOIZED) 

0.6 

KEV/CHAN 

CNEL  ' OOJu  SEC.  / 

CATA  ITITLE(1,P1) /b)HFTG 

• 

RA 

3EF 

n 

21 

(JNOXIOIZEO) 

0.6 

KEV/CHAN 

CNEL  ‘ OC  3G  SFC.  f 

CATA  1TITLE(1,E2)/C>0HFTG 

• 

AM 

241 

ORT 

EC  1/21/69 

0.6 

KEV/CHAN 

CNEL  1CG0C  SEC.  / 

DATA  ITITLE(  l,r  3) /blHFTG 

• 

AM 

241 

ORTEC  1/21/69 

0.6 

KEV/CHAN 

. CNEL  10U0C  SEC.  ' 

DATA  ITI7  LEU,  E-> /C»  JHFTG 

• 

AM 

241 

or  r 

EC  1/21/69 

0.6 

KEV/CHAN 

CNEL  10C0G  SEr.  / 

DATA  ITITLEU.^5) /E  1HFTG 

• 

AM 

241 

ORTEC  1/21/69 

0.6 

KEV/CHAN 

CNEL  3CCJC  SEC.  / 

DATA  ITITLE(l,rG) /feOHFTG 

• 

RA 

3rr 

it 

15 

(JNOXIOIZEO) 

UNBIASED  RUN 

C *.  3 , 0 C 3 SEC.  / 

DATA  ITIT LE( 1,67) /6THFIG 

• 

PA 

3 E f 

it 

15 

(JNOXIOIZEO) 

UNBIASED  RUN 

C AO, OLD  SEC.  / 

DA  I A ITITLE(l,F3)/bJHFTG 

• 

RA 

3 E r 

it 

15 

( JNOXIOI7ED) 

0.6 

KEV/CHAN 

CNEL  * 0 00 C SEC.  / 

DATA  I T I T L E ( 1 , D 1 ) / 6 T M F I G 

• 

RA 

3 E f 

it 

15 

( JNOXIOIZEO) 

0.6 

KEV/CHAN 

CNEL  l GOOD  SrC.  / 

DATA  ITITLEd,  bO» /63HFTG 

• 

RA 

3Er 

n 

14 

(JNOXIOIZEO) 

0.6 

KEV/CHAN 

CNEL  ‘ COTl  SEC.  / 

DATA  ITITLE(l,bl) /6THF-G 

• 

RA 

3Er 

it 

1* 

( JNOXIOI7EO) 

0.6 

KEV/CHAN 

CNEL  ‘0,3-0  SF3./ 

DATA  ITITLE(l,o?) /GTHFIG 

• 

RA 

3EF 

n 

13 

( JNOXTniZEO) 

0.6 

KEV/CHAN 

CNEL  * C , 3 J 0 SEC./ 

DATA  ITT! LE(l,b3) /GO^FTG 

• 

RA 

3EF 

n 

13 

( JNOXTOT7ED) 

0.6 

KEV/CHAN 

CNEL  ‘ 0,100  SEC  ./ 


DATA  riTL-(l,0^!/6HFIG 
CNEL  * C ,3  JO  SEC  . ' 

. RA 

9 i r 

if  12 

( )NQXTDI7ED) 

0.6  KEV/CriAN 

DATA  ITITLEd,  6r>) /bOHrjr; 

CNEl  i C ,0c  0 SEC  . / 

. PA 

DTP 

ft  17 

( JN0XI0I7ED) 

O.b  KEV/CHAN 

TATA  ITITLE(1,S'>) /50HFIG 
CNEL  d,  I'JC  SEC./ 

. RA 

3Er 

If  11 

( JNCXIOI ZED) 

0.6  KEV /CHAN 

DATA  lTITLEd»c7)/53HFIG 
CNEL  ‘ G,  uC  C SEC  . / 

. PA 

3 E r 

If  11 

( JNOXIDIZED) 

0.6  KEV/CHAN 

DATA  ITXTLF(  i,S1) /GOHFTG 
CNEl  '0,300  SEC. / 

. RA 

3E-“ 

If  13 

( JNOXIQIZEO) 

0.6  KEV/CHAN 

DATA  ITITLEd,  b3> /&0HFTG 
CNEL  * G,  ..0  0 SEC . / 

. RA 

3EP 

if  13 

( JNOXTOIZED) 

O.b  KEV/CHAN 

DATA  ITITLE(1,70) / 60  HFTG 
CNEL  j 0,3;  O SEC.' 

. AM 

2 + 1 

CRT 

EC  i/d/6y 

O.b  KEV/CHAN 

DATA  I^ITLEd.Tt) /S3HFIG 
CNEL  10 ,0L  3 S EC . / 

. AM 

2d 

ORTEC  1/21/69 

O.b  KEV/CHAN 

DATA  1TITLE(1,72) /61HFTG 
CNEL  It,  JOG  SEC./ 

. P'J 

290ft2  A ? 1 1 N I 77 

0.6  KEV/CHAN 

DAI  A 1TIT LE( 1,73) /60HFIG 
CNEl  1 L , 3 u C SEC./ 

. PU 

29  0 

if  2 

ARMANI  / 7 

0.6  KEV/CHAN 

CAT  A I7ITLEd,74l /&OHFIG 

cnel  rc,ooo  ?fo./ 

. PM 

29  J 

« 1 

SRMA-I  77 

0.6  KPV/CHAN 

DA  1 A lTITLEd,  75) /oOHFTG 
CNEL  ? t: , 3 0 C SEC./ 

. PU 

290 

* 1 

ARMANI  7 7 

O.b  KEV/CHAN 

DATA  ITITLE(i,73'/el-iF7G 
CEL  tr,0Gu  SEC.  / 

RA  D 

E r <f 

10 

( 3XIDI ZED) 

0.6  KEV/CHANN 

DATA  lril LE(J,77) /53MFIG 
CEL  N C , 0 30  SEC.  / 

RA  1 

Er  n 

10 

(OXIDIZED) 

0.5  KEV/CHANN 

DATA  ITITLE(1,7A) /63AFTG 
CEL  i 0,  COO  SEC.  / 

PA  D 

E-  if 

11 

( 3XIOI ZED) 

0.6  KEV/CHANN 

DATA  1TITLE (1,79) /50HFTG 
CFL  AC,  CO  3 SEC.  / 

RA  D 

E 7 * 

11 

(OXIDIZED) 

0.6  KEV/CHANN 

DATA  I7ITLE(1,8(?)  / o 0 H c I G 
i CEL  < 0,030  SEC.  / 

PA  D 

E * if 

12 

(OXIDI ZED 

0.6  KEV/CHANN 

DATA  ITIT  LEd,  31) /60HFTG 
CEL  A C , 0 o r.  SEC.  / 

RA  n 

E-  ft 

12 

(OXIDIZED) 

0.6  KEV/CHANN 

OAT  A ITITLE(1»32)  / oO  -l  CIG 
CCL  **C,  0 0 0 SEC.  ' 

FA  7 

-r  * 

13 

(OXIDI 7ED ) 

0.6  KEV/CHANN 

DATA  ITITLEd,  87) /&0HFTG 
CFL  AC,  03  3 SEC.  / 

RA  D 

E - ft 

13 

(OXIDIZED) 

0.6  KEV/CHANN 

DATA  ITITLE(l,o4) /SOHFJG 
CEL  l f , 0 3 C SEC.  / 

PA  7 

E r >f 

14 

(OXIDI ZED) 

0.6  KEV/CHANN 

DATA  ITITLE(1,6")/S0HFTG 

cfl  oc,oo:  SEC.  / 

RA  7 

" ft 

1A 

(OXIDIZED) 

0.6  KEV/CHANN 

DATA  ITITLE(l,a<i) /SOHFEG 
CEL  * C, COl  SEC.  / 

RA  7 

E - ft 

15 

( OXIDT 7ED) 

0.6  KEV/CHANN 

DATA  I7ITLE(i,67) /50HFIG 
CEL  Hi, COO  SEC.  ' 

RA  0 

E it 

15 

( OXIDIZFD) 

0.6  KEV/CHAnN 

nftTA  IT27LE(1,83>/5(HF7G 
CEL  h c , C 0 0 S£D . / 

RA  7 

E-  ft 

lb 

(OXIDI 7ED) 

C.6  KEV/CHA9N 

DATA  ITI iLE(l,  39) /50HFTG 

C CL  AC,jOO  SEC.  ' 

RA  D 

" ft 

16 

( OXIDISED) 

0.5  KEV/CHANN 

DATA  ITITLrd,  90) 

CHANNEL  10  K5  EC  . / 

. AM  2 

91  (ORTEC 

) UHqT  A SEO  RUN  0.6  KEV'C 

DATA  lTITLEd,  Fl) /GQHFTG 
CHANNEL  1C  <S  EC  . / 

. AM  ? 

♦ L (ORTEC 

) UNBIASED  RUN  0.6  KEV/C 

OATA  ITITLE(1,92) /609FTG 

' . 9 A1' 

74 

<5  ROUND  C 

OJNT  ( °I AS  ED ) hOOOO  SEC 

DO  O O O O O O 


c.  / 

DATA  I7ITLE(1,93)  /i)HFTG 
C.  ' 


BAG  < j R3  UNO  CDJNT  BIASED)  4 0000  SEC 


,»»***♦»*».»»»  ♦ ■.***** 


C the  PIFST  S°ECTPI.'1  TO  BE  PIOTT-D  MUST  BE  RES  I 5 T EFED  BY  THE  USE  OF  THE 
C SKIPF  COUPOL  FARO 

C NTAG  (1  ,N)  IS  THE  TAGWO-O  TO  3E  P.DTTEO 

C NTAG(4,N)  IS  the  LAST  CHANNEL  TO  BE  PLOTTED 

C NTA  GC3,\)  IS  "HE  rIRST  CHANNEL  T j P£  PLOTTED 

C NTAG  (2, N)  IB  T UE  1‘VPE  OF  = L3T:  1 = St-ALL  FOR  PUBLICATION 

C 2 = LARGE 


C ♦-  '* 


*♦*-**• 


»»••»»*  + » » » ► » » * • »»**♦*.* 


READ-  ,YA(l),YA(l),vA(?),XA(2) 

CIRST,  THE  TAG^mns  T0  BE  SLOTTED  AND  OTHER  DATA  ( NTAG ) 
ARE  FEAo 


J=  u 

5 J=J  + 1 

pEAO'  , (NT AG ( I, J) ,1=1,4  ) 

IE (NT  AG  ( 1 . J)  .GT . f ) GO  TO  5 
J = J-1 

PRINT-’,"  TAGW3RDS  ”,  ( NTA  G (1 , JJ)  , J J=  1 , J ) , “ TO  BE  FLOTTED 
TALL  AkYODR (NTAG, J) 

N=  J 

DO  9 C C I = f,N 
YAXIS^XA (NTAG( 2,1) ) 

YAXIS  = YA ( NTA  G ( 2, I ) ) 

A PL  E = XAXIS  *■  4.0 
K = NTA  G(«t,I)-NTAG(3,T)  +1 
BEGIN  NOFXING  ON  SPc:TRA 
11  NN=1 
K’K  = 1 

>-  ^•*■4*.  *-****♦»  ♦*«•*■»■>  * + +4  + +*  +» 

READ  IN  THE  X VECTOR 
40  BUFFER  I J ( 1 , 0 ) < HW ( l)  , M « ( lb C ) ) 

KK  = lJf IT ( 1 ) 

CALi.  EC'K$ET  COUNT,  IT  O) 

5 0 DECOrE(l5C,f>«5,MW)  M,  HH,7 

9 0 FGrNAT  (2X  ,1-4 ,2X,  A - ,?.TFf  . 0,9  (/25rS.  0)  , /,3r5.  0) 

PRINT  * t“  taghur:  NIHPFR  ",M,"  BEAD 
IF (1  -NTAG (1,1)1  11,70,903 

70  CONTINUE 

IF (KCUUT.GT . 0)  GO  T1  1°0 

71  CONTINUE 

DO  ICC  J=l, 2Bn 
Y ( NN ) = 7 ( J ) 

10  0 NN=NL ♦ 1 

IF(<f.nE.0)GO  TO  40 


SET  UP  THE  X VECTOR 

X ( 1)  = cl0AT  <NTA0<  3,1 ) ) 

no  10  1 1 = 2 , < 


io  xnn  = x(ii-i)M.r 

x(kh  i ) = X(l) 

X ( K 2 ) = float  (<) /XAXTF 


O O O O O O O O i *■>  o o 


4»*»*4*'  »•»*»•*•  *.*•»♦.*♦* 

EppQF  EDIT  3 EE  PROGRAM  33RT°T 
IP  IE  EOl'CLS  7EPO  THAN  THE^E  acE  ERRORS  IN  HE  V VECTOR  FROM  SCRIPT 
NE  IS  THr  K"JM PER  ERRORS  IN  THE  VECTOR  Y 
FEAD  1 C 3 0 , 1 E , N E 
Ic  (IE.  ED. ti' CO  TO  3 0 G 

*■***■♦*»*»■•*♦*.♦*****■*  *♦♦•*#♦*«.».»*****♦*»•.»*».  ****»» 

SET  UF  PLOTTER 

CALL  PLOTtC. 0,-3. 0,-3) 

IF (NT AG( 2,T)  .GT  .1 ) TO  TO  20 
''all  plot  ( i.  o,  7,r , •»' 

CALL  PLOT  (1.0,6. r.,  2> 

CALL  PLOT  (<>.0,9.;  , 3) 

CALl  PLOT ( 9 • 0 » Q » • , 2) 

CALL  PLOT (12.0, 3. C, 3) 

CALL  PLOT  (12.0,3.0,?) 

CALL  cl. Of  (9.  u , 1 . C , 3) 

CALL  FLOY U. 0, l.r , 2) 

CALL  PLOT  (2.23,2.0,3) 

CALL  PLOT <2.23, C. 3,?) 

CALL  PLOT (11. 0,3. 3, ,) 

CALL  PLOT (11.3, 2.0, ?) 

CALL  PLOT (2.23 , 2.C ,?) 

CALL  PLOT (2. 73, 2. 5, -3) 

CO  TC  21 

190  PRINT’,"  PAD  OATA  READ,  TASWORO  * m 
CO  TC  71 

20  CALL  FLOT ( l.C, 2.0, -3) 

21  CONTINUE 

SCALE  Y VECTOR 
2 3 C OC  21<?  J=1,K 
210  Y ( J ) = Y ( J + NTA G (3,1) +N OFF-1 ) 

T=P=Y( 1) 

DO  2 2 L J = 1 . K 

IF ( Y ( J) ,GT.T)T=Y(J)  x 

I F ( Y (J)  ,LT  .3) 3 = Y ( J ) 

220  CONTINUE 
Y (K+3 ) =9 

Y ( K+  2 ) =(T-«)  /YAX1S 
PLOT  PFS'JLTS 

CALL  4XIS(0.u,0.0, 7HCH ANNEL , -7,YAYI3,0.G,Y(K+l),X(K+2)) 

CALL  A X I S ( 0 • 0,  fl.O,  NH3DMNTS,6,  YAXI  5,  9 0 . 0,  Y ( < *■  1 ) , Y ( K +2)  ) 

CALL  LINE  <X,Y,  <,1,  [!,  1) 

PPINT  1033, M 

CALL  SYMBOL ( 1. U,*YM, 0 .10F,ITITLE< 1 ,M>  ,0.0,30 
CALL  FLOT  (ABLE, 3. 0,-3) 

300  CONTINUE 

CALL  PLOTc 
STOP 

300  PO  31L  J=1,NE 

PE*P  1C1P  , IC » AC 
I C = I c *1 
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TIP  V ( IC ) = AC 
GO  TO  2H 
10C0  FOf-  < 2 1 E> ) 

1010  FORUM  (I'.‘,F10.0) 

1020  FOt MCT(2ri0.0) 

10  30  FORMAT (IX ,"TAGWO-n  ">L  0"  T ED” , T 18 , [ 5 ) 
FND 


GUPPCUTIME  AF  YCP~  (MARFA  Y,N) 

OlnEMSI  ON  MAR? AY  ( + , >1 ) , T S A V(4  ) 

MF=N-1 

00  Of  I = 1 , N P 
t'sIM 

00  5T  J = K,N 

Tc (MURRAY  (i, I)  ,LE.MA°FCY ( 1 , J ) ) GO  TO  50 

00  30  m=i,k 

IGA'/(*)  - HAR?AY(m,T) 

NARRA Y (M, I)  = NAR?AY(M,J) 

30  M ARRAY (Y » J)  = ISA7(  <) 

50  CONTINUE 

f,‘  F T U F u 

END 


?(JFF  CUT!  ME  HALVE(SP-C,  rWHM,  RHHrO 
REAL  LWHM,RWHM 
DIhENSIUN  SPEC(2100) 
r<RAX=lJG  . 

DO  10  0 1 = 10,2100 

IF(SPECd)  .Lr.PKYAY)  GO  TO  100 

DKA  A Y = S'3  EC  (I) 

IPEAYrl 

ICO  CONTINUE 

HA  FMA  X=p!<, MAX/2  • 

, 0 = 1°  F A K 

lr0  J=J+1 

2 C 0 IF  (SPEC  ( J)  .GT.HA-MAXO  co  TO  150 

9WH-  = J + (ROEC  <J)  -Mfl  r-iA  y)  / (S»EC  ( J-ll  - 3P.EC<  J)  > 

J= I p E A Y 
2r  C J = J-1 

IF  ( $ P E P ( J ) .GT.uArMAX)  '■-0  TO  250 

LWHM=j  + (HA^NAX-GFEC  (.)))/  <S3EC  (Jn  > - SPEC(  J)  > 

rwwr)  = pwhv-LMHM 

PRINT*,"  PK*AX  = ",'>KMAX>"  L WH  M = FWHM  = ",FWHM, 

C,RWUf 

PETMFN 

FND 
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RWHM=  •’ 
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Appendix  E 


User  Instructions  for  SPECTRE 

Program  SPECTRE  is  designed  to  allow  the  user  to  select  the 
tagword  to  be  plotted,  choose  a large  or  small  format,  and  select  the 
range  of  channels  for  the  plot  with  one  card  for  each  spectrum  to  be 
plotted.  A set  of  title  statements  are  initialized  by  the  use  of  DATA 
statements . 

The  file  containing  the  spectra  to  be  read  must  be  placed  on  file 
TAPEl  as  for  Program  SCRIPT,  but  the  first  spectrum  need  not  be 
registered.  It  is  more  efficient  to  advance  the  file  pointer  to  just 
before  the  smallest  tagword  to  be  plotted,  but  this  is  not  absolutely 
necessary. 

A series  of  2N+1  data  cards  are  required  to  plot  N spectra.  The 
first  card  must  have  the  axial  dimensions  in  inches  for  the  two  plot 
formats  in  the  following  order:  counts  (small),  channel  (small), 
counts  (large),  channel  (large).  These  are  normally  5.0,  7.0,  8.0, 
12.0.  In  the  large  format,  the  counts  axis  is  limited  by  the  width  of 
the  plotter  to  less  than  nine  inches  but  the  channel  axis  is  unlimited. 
The  small  format  dimensions  should  not  be  changed. 

Next  should  come  a series  of  cards  with  the  information  for  each 
spectrum  to  be  plotted.  Each  card  should  contain  in  order  the  tagword 
to  be  plotted,  the  type  of  plot  (1  = small,  2 = large),  and  the  first  and 


last  channels  to  be  plotted. 

Finally,  error  data  cards  are  added  for  each  spectrum  to  be 
plotted.  Each  card  should  contain  IE  and  NE  in  215  format  with 

IE  = Is  there  an  error? 

0 = Yes 
Not  0 = No 

NE  = Number  of  errors 

If  IE  = 0 there  should  be  NE  cards  containing  the  channel  number 
and  new  value  for  each  error  in  format  15,  F 10.0.  These  will  be  read 
before  the  next  error  data  card. 

The  order  of  the  tagword  cards  is  unimportant  since  they  will  be 
placed  in  sequence  from  smallest  to  largest  by  subroutine  ARYODR. 
However  the  error  edit  cards  must  be  in  correct  sequence  to  match 
the  tagwords  after  the  tagwords  are  ordered  by  ARYODR. 

The  SCOPE  control  cards  required  are  similar  to  those  for 


SCRIPT. 


Appendix  F 


Overview  of  Program  ALFAIC 

I 

This  appendix  is  a discussion  of  the  general  flow  of 
action  in  program  ALFAIC.  Refer  to  the  block  diagram  in 
Figure  24. 

The  main  program  first  reads  certain  controlling  data 
including  the  run  TITLE  card,  I*AMELIST  CONTRL,  and  the  back- 
ground values  BACK,  which  are  then  sent  to  ARYODR  to  be 
placed  in  correct  sequence.  The  program  then  calls  RDSPEC 
to  read  in  a reference  spectrum  from  file  TAPEl.  Subroutine 
HALVE  computes  the  FWHM  for  this  spectrum.  The  program  then 
calls  SPCTRM,  which  in  turn  calls  RDSPEC  to  read  the  spectrum 
to  be  analyzed.  Then  the  tentative  positions  of  peaks  to  be 
resolved  are  read  into  PEKPOS  and  sent  to  ARYODR  for  ordering. 
On  subsequent  iterations,  the  program  will  read  in  a speci- 
fied number  of  reference  spectra  to  use  to  analyze  this  same 
unknown. 

Now  the  program  calls  BGRND,  which  can  perform  one  of 
three  options.  BGRND  can  automatically  compute  the  back- 
ground using  subroutines  LSQPOL  and  MATINV.  It  can  also 
compute  an  interpolated  background  function  through  a 
designated  set  of  points.  The  third  option,  which  was  used 
throughout  this  study,  is  to  compute  a number  of  regions  of 

80 


REFPEK 
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Fig.  24.  ALFAIC  Overall  Flow  Chart 
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constant  background  from  data  in  BACK. 

Next,  the  program  calls  REFPEK  which  computes  certain 
parameters  for  the  reference  peak.  These  may  instead  be 
specified  and  the  program  will  compute  a reference  to  match 
or  the  program  can  compute  a reference  from  the  spectrum  to 
be  analyzed.  The  parameters  are  based  on  the  highest  and 
lowest  nonzero  values  in  the  reference  spectrum,  the  maxi- 
mum value  in  the  reference  peak  and  the  one -third  maximum 
point  to  the  right  of  the  peak. 

The  program  then  checks  the  peaks  to  determine  if  they 
can  be  fitted  within  the  range  of  the  spectrum  using  a 
reference  peak  of  the  length  computed  in  REFPEK. 

Finally,  the  program  calls  RESLTS,  which  passes  the 
data  on  to  PREPAR.  When  the  computed  spectrum  is  returned 
from  PREPAR,  RESLTS  computes  the  areas  and  estimated  errors 
for  each  peak  and  the  total  areas  of  both  the  unknown  and  the 
resolved  or  calculate  spectra  and  the  square  root  of  the 
sum  of  the  square  of  the  difference  between  these  two  spectra 
at  each  point. 

PREPAR  sets  up  the  problem  for  DAVIDN,  which  calls  FCN 
and  then,  as  required,  calls  READY,  AIM,  FIRE,  DRESS,  ORDER, 
and  FCN  again.  FCN  utilizes  SHIFT,  FUN,  MATINV,  and  DERIVP 
to  help  it  search  for  the  minimum  Chi -Squared  Function.  These 
routines  collectively  do  the  work  of  calculating  a spectrum 
which  is  the  result  of  a constant  for  each  peak  multiplied 
by  the  reference  spectrum  at  each  point  relative  to  that 

peak  location.  These  individual  components  are  then  added  to 
give  the  combined  spectrum. 


Appendix  G 


Input  Cards  for  ALFA  1C 

Introduction 

-ALFA1C  allows  the  user  a plethora  of  options,  but  placing  the 
input  cards  in  the  correct  order  is  critical  if  the  program  is  to  func- 
tion as  desired.  The  program  can  be  run  in  two  main  ways,  analyzing 
real  data  or  internally  generated  data. 

Conjur  Option  Conversion 

As  listed,  the  program  analyzes  real  data  from  file  TAPE1  and 
the  routines  CONJUR,  FACTOR,  PPOISS,  PGAUSS,  and  OUTPUT  are 
not  used.  By  replacing  the  calls  to  RDSPEC  with  calls  to  CONJUR, 
the  program  may  be  converted  to  analyzing  randomized  Gaussian 
spectra.  Input  6 should  be  substituted  for  Input  5 if  this  version  is 
used. 

Standard  Version 

The  unknown  spectrum  and  NRSPEC  reference  spectra  should  be 
placed  in  file  TAPE1  in  the  following  order: 

Reference  spectrum  1 
Unknown  spectrum 
Reference  spectra  2-  NRSPEC 

FILE  and  LDSET  control  cards  arc  needed  as  for  Programs  SCRIPT 
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i 


r 


I 


and  SPECTRE. 


Input  Cards 

Input  1.  Title  Card 
FORMAT  (8A10) 


The  title  for  this  complete 
sequence  is  read  from  the  first 
card.  End  of  file  or  blank  termi- 
nates the  run. 


Input  2.  NAMELIST  CONTRL 
NBACK  = 0 

= +n 

= -n 

Note  1:  im  < 50 


Background  determined  auto- 
matically by  program. 

Background  option  I.  The  X and 
Y coordinates  for  n background 
points  will  be  read  in  i program 
calculates  a curve  passing 
through  all  points. 

Background  option  II.  n back- 
ground intervals  will  be  read  in 
(background  is  a constant  between 
the  limits  of  each  interval). 

If  n > 50,  the  program  will 
truncate  it  to  50. 


Note  2:  n 5 3 for  background  Option  I (positive  sigii). 


NPKRD  > +1 


Initial  peak  positions  to  be  read, 
as  values  of  channel  number. 


NREF  = 0 


= -1 

NSSW1  = 1 


Reference  peak  to  be  selected 
from  the  spectrum.  (This  option 
is  generally  not  recommended.  ) 

Reference  peak  spectrum  to  be 
read  from  tape. 

Print  the  results  from  variable- 
metric -minimization  routines. 


= 0 


Suppress  the  above. 
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AD-A064  056 


UNCLASSIFIFD 


AIR  FORCE  INST  OF  TECH  WR I 6HT -PATTERSON  AFB  OHIO  SCH— ETC  F/6  7/2 
COMPUTER  CODE  TO  ANALYZE  ALPHA  SPECTRA  USIN6  A SPECTRAL  STRIPPI— ETC(U) 
DEC  78  J R HARSTINE 

AFTT/GNF/PH/78D-16  NL 


1 


MICROCOPY  RESOLD  I ION  II  SI  CHARI 


I ( 

( v 


NSSW2  = 1 


MAXREF  = 


KEPREF  = 0 


KEPDAT  = 0 


STDEV = 


List  the  raw  data,  resolved  peaks, 
and  composite  spectra. 

Suppress  the  above. 

The  maximum  number  of  channels 
in  the  reference  peaks. 

Read  in  NRSPEC  new  reference 
spectra. 

Use  the  reference  spectra  from 
the  last  analysis. 

Read  in  a new  spectrum  to  be 
analyzed. 

Analyze  the  last  data  spectrum 
again. 

The  standard  deviation  of  the 
Gaussian  spectra  to  be  analyzed. 


STDEL  = The  amount  the  different  refer- 

ence spectra  are  varied  around 
STDEV. 

NRSPEC  * The  number  of  reference  spectra 

to  be  used. 

it  3.  Background  Data  Cards  (not  used  for  NBACK  **()). 
DB(I ) , BACK(I) , 1*1,  NBACK 

FORMAT  (16F5.0) 


DB(I)  = 


11-15,  etc. 


For  NBACK  > 0 (option  I)- 


DB(I)  is  the  channel  number  at 
which  a background  value  is  to 
be  specified. 

For  NBACK  < 0 (option  II): 

DB(I)  is  the  channel  number  at  the 
lower  end  of  the  Ith  interval  for 
which  a background  value  is  to  be 
specified.  DB(1)  is  supplied  by 
the  program  and  may  be  left  blank. 
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6-10  BACK(I)  = The  value  of  the  background  for 

the  Ith  background  position 
(NBACK  >0),  or  for  the  1th  back- 
ground interval  (NBACK  <0). 

Note  1:  The  pairs  DB(J),  BACK(I)  may  be  given  in  any  order.  The 
program  will  rearrange  them  in  increasing  order  of  DB(I). 

Note  2:  For  NBAC1O0,  the  program  requires  backgrounds  for  the 
first  and  last  channels  of  data.  If  not  given  explicitly,  the 
program  will  use  BACK(l)  for  the  first  channel,  and  BACK 
(NBACK)  for  the  last  channel.  These  program-supplied  values 
may  be  counted  for  meeting  the  minimum  requirement  of  3 
background  points. 

Input  4.  Peak  extraction  channel 

IPK  = The  desired  channel  number  of 

the  maximum  count  in  the  512 
channel  segment  to  be  analyzed. 


Input  5.  DO  Loop  Sequential  data 

Five  sequences  of  these  cards  will  be  required  for  each 

analysis. 

Input  5a.  Reference  Peak  Parameter  Card  (Use  only  if  NREP=0. ) 
STOIPR,  DELTA  1 , DELTA2,  DELTA3 


FORMAT  (4F10.  5) 

I- 10  STOIPR  = 

II- 20  DELTA  1 = 

21-30  DELTA  2 = 

31-40  DELTA3  = 


The  channel  number  of  the  maxi- 
mum count  in  the  reference  peak. 

Channels  from  the  front  (high- 
energy)  edge  of  the  reference  peak 
to  the  "third-height"  position  on 
the  front  edge. 

Channels  from  the  back  (low- 
energy)  edge  of  the  reference  peak 
to  the  third-height  position. 

Channels  from  STOIPR  to  the 
third-height  position. 
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Note:  For  MCA  operation,  the  "third -height"  position  is  defined  as 
STOIPR  and  hence  the  program  sets  DELTA3=0.  0. 


Input  5b.  Peak  Position  Cards  (Use  with  NPKRD  > +1. ) 

PEKPOS(I) , IFIXPK(I) 

UNFORMATTED  (LIST  DIRECTED) 

PEKPOS(I)  = The  estimated  third-height  loca- 

tion of  the  Ith  peak  to  be  considered 
by  the  program,  specified  in  chan- 
nel number. 

IFIXPK(I)  - 0 or  blank  The  program  will  adjust  PEKPOS(I) 

for  the  best  fit  to  the  data.  (Stand- 
ard option. ) 


> 1 


PE'KPOS(I)  will  be  held  fixed  by 
the  program. 


Note:  The  cards  may  be  in  any  order.  The  program  will  rearrange 
them  into  ascending  order  of  PEKPOS(I), 

Input  5c.  Peak  Position  Termination  Card 

PEKPOS(l)  < -1.0.  This  will  terminate  peak  locations  and 
the  program  will  process  a new  set  of  data,  beginning  with 
Card  1. 

Input  6.  Conjured  Data  Sequence  Cards 

These  are  used  instead  of  INPUT5  when  running  the  program 
using  data  to  be  created  by  Subroutine  CONJUR  and  the  associated 
functions.  The  two  statements  in  ALFA  1C  and  SPCTRM  calling 
RDSPEC  must  be  replaced  with  calls  to  CONJUR  for  these  options. 

Input  6a.  Channels  of  desired  peak(s) 

IPEAK(N) , N = 1 , NPEAK 


UNFORMATTED  (LIST  DIRECTED) 


IPEAK  = 


The  center  channel  of  the  Gaussian 
peak  to  be  created  by  calling  func- 
tion PGAUSS. 


Input  6b.  Heights  of  desired  peaks 
IHITE(N),  N = l,  NPEAK 
UNFORMATTED  (LIST  DIRECTED) 

IHITE(N)  = The  approximate  height  in  counts 

of  the  peak  in  channel  IPEAK(N). 

Input  6c.  Peak  position  card3 

These  are  the  same  as  Inputs  5b  and  5c. 
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PROGRAM  ALFaiCd1)0  JT  ,ruTP"T  ,TA3: 5 = [ NPUT , T 1 P E6  =0UTPMT  , PUNCH, 

1 TAFE7= ^UNCH, TAPER. rArEl) 

F'-tG*  A * AlP^A-TT  HAG  OrTAlNEO  F?DW  PROGRAM  AUTPFIT  WRITTEN  RY 
J.F.  COMFORT  IF  f RGOMNr  NATIONAL  -AlORATORY 
THIS  PROGRAM  ALLOTS  T Hr  FOLLOWT‘J-,J 
A MAXIMUM  Op  21:11  CHANNELS  FOR  A R E"E  R ENG  E PCAK 

A MA  y I M J M OF  2101  GUNNELS  FOR  A >PECTRJ* 

A rA  Y1MJN  OF  ;d  GMAJNPLS  for  a nGKSROUNO  TNPJT 

COMMON  /M  1(  = 42)  ,*SSW1,NSSW? 

COmMCN  7 A 0 2/  lUHNSdRi  ) , STOP'/ 1 ( ? 1 0 0 ) , ?SWX(210C>,  ITIMES,  NNX 
COMMl  • TITLE! 21)  , FEKPOS  (?0),rjrrXP'<(20>  , BET  A , GST  ATE, 

1 GSF\/Al,  NT,  NDA.C,  KS,  LO 

COMMON  /;  0,/  INT;  rF  ( , ltj ) , ST  ORrr  C 2 1 P 0),0F.Tfil,DELTA2,O£LTA3, 

1 STC1PR,  REF  1 A X , 111,  1 02,  IDT,  'll,  NE»  WREF 
Co  MM  ( N /A  0r»  / lAC^2(210d  , 0*0(50),  3ACK<EC)»  NJ,  N3ACK,  NPKRO 
COMMON  / 0 A T A / ni5T(?10r),  0 0 UNT  T ( 21 3 0 ) , RJWIO,  TODAY,  FRED,  OFTEN, 
1 r. OFTEN,  NCHANI,  It,  TJ,  ISPECT,  I-MT,  IER R , M« XF EF, MAXP2 
COMMON  /CiPTN/  <E  = RE-  , KrPOAT 
nihFf  SIUN  REFOTOlEn-*,  10)  ,cWMM(m?> 

FOUI  VALt.'CE  (IRLNK,  A*.  N'/,  DDL  NK) 

"*AT  A ( GnLNF=6u  ) 

NAMEL1ST/CONTRL/NPAO< ,HPKRD,NPEr, NSSW1,\|>SW2,MAXREF,KEPREF,KEPDAT 
C,3TnEV,STO£L,NPSFFC 

...START  CALDOMP  AND  GET  OAT"... 

1EKF  = 0 

CALI  OAIE(TCOAY) 

GO  U DO 

. . .TERMINATE  PROGRAM... 

4R  WPITF  (t,  101) 

101  FOkM/T  (3X.18H  Ef:0°L0T  EXECUTED) 

STOP 

...READ  ANO  WilTR  TTJLF  CAPO.  EXIT  ON  E.O.F... 

50  PEA0(‘,1P)  (TITLE(I) , 1=1, 3> 

I F ( F OF ( 5 ) ) *H,43 

49  TF  (TITLE(l)  .EO.r-.N<)  GO  TO  L5 

WPiTF(G,30)  ( T IT  L =■  (T>  ,1=1,6  )» TODAY 

. . .aNFIJT  0"  DATA... 

FFMn  <C  ,CCNTP'_) 

IF ( EOF (5 ) ) 1*5 ,51 

51  PRINT  CONTRL 

yEPKFF  IS  AN  OPTION  TO  ALLOW  THF  JSER  TO  <EEP  THE  SAME  REFERENCE 
PEAK  FCR  ALL  AnA.YSIS 

•YEPREF.GT.O  ALLOWS  THE  USE  OF  THE  SAME  REFERENCE  PEAK  FOR  ENTIRE 
ANALYSIS 

KEPaFF.ED.J  -EO')IRrS  A NEW  reperenge  PEA<  WITH  each  spectrum 
ANAL V7E0 

Kc-f-J*  T *>  IN  OPTION  TO  ALLOW  TH-  SA1E  DATS  SET  TO  RE  USEO  FOR  ALL 
ANALYSIS. 

tr-rppeT  • GY « 0 ALLOWS  "HF  SAMF  DA”  A SET  TO  BE  USED  FOR  THE  ANALYSIS. 
*EPDAT  .EO.O  REDJIRE5  A NEW  DATA  SET  FOR  EACH  ANALYSIS. 

M A XP  2 = NAXREF  ♦ 2 
IF  (NPPF.GT.il)  VRCP  = 10 
IF  ( Na  ACK ) TG,11(,A1 
70  M J = -NRACK 

Na AC K t 2 
GO  TO  90 
«0  NJ  = NrACK 


» 
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fj0£0  = 1 

<Q0  TF  (NJ.GT.bO)  NJ  = on 

REAP <5, 14)  (Cn(I),'UC‘MI>l!=l,NI> 
no  *1  1= 1 « nj 

ni  print*  ,*•  m cm,”)  - ",n9(i)f"  = mackci) 

IF  U'J.ST.l)  CAL.  APYOGR<Or,RAC<f YJ,  2) 

C . . .INITIAL! 7ATICM.  . . 

inn  CONTINUE 
LOrQ 
NAOJ=C 

CT  01 FR  = 0.3 
on  115  1=1, 20 
i cb  iFixPKfn  = n 

DO  110  I=l,MA<PEr 
110  STOrFFO  = 0.0 
C . . .MO*  E jAT a IMP  IT  . . . 

T c I f K CF ) i: j, 130, 1J0 
1?0  IP ( KFFR  £F • GT .3)00  TO  7 '8 

PRI N’*,”  CALLING  PP>OFC  POP  REEMN3E  PE  V < " 

00  ltC  l I J<=l,Nr?JPE: 

fall  kPSpECISTOPEF  , N * NT  AG»“AXR£CI 

K-c-jot  = u 

00  77  I = 1 » N 

PEFSTCC,  I JO=STOvE- (I) 

77  CONTINUE 

call  H Al  V E ( R ECSTD ( 1 * X J ^ ) » FWH M ( I I < , 1 ) ) 

FWHM  (IK  » 2)  = IK 
1001  CONTINUE 

GO  10  lJ'J 

▼7  8 DO  95  I = 1 » NS  T 0 R 

STOP  E F (t) =PEF5T0(I ,1 JK) 

99  CONTINUE 
100C  CONTINUE 
GO  TO  15  0 

13F  PEAEM5»l.i>  S TOI°? , OIL  TA  1 , DE L TA2 , 0 El  T A 3 
IE  G CONTINUE 

C ...GET  A ST  AMD A RD  SPECTRUM  DATA  SET... 

JF (X  FLOAT  .GT  .0 ) GO  TO  175 
170  FALL  5PCTRM 
175  IF  (JERR.GT.O)  GO  TO  L ' 

WRIT  F (f  , 7i* ) NOHANN 
CELT  A3  =0.0 

C ...READ  PEAK  POSITIONS,  DISTMOE 

200  DC  220  1=1,21 
210  REAO-  ,RE<POS(I)  , lrIV»^(I) 

IF  trfKPJSm  .LE.-l.  ) GO  TO  230 
PRUT",”  PEAK*  ”,I,”  = ",cEKPOS(n 
220  CONTINUE 
NT  = ?r 
GO  TF  24 G 
230  fK  = I - l 

20  IF  (f.T.&T.l)  CAL.  ARYOr’R(PEKPOS,IrI<PK,MT  ,2) 

I c Ki  = PEKr  OS ( 1 ) 

330  00  3!  0 I = 1,NCHANN 
S^c  F-AWYd)  - CO'JMTS(t) 

C ...CALC.  dACKGRO’IMD,  SELECT  RErERENC£  3E  A K AND  PICK  PEAKS... 

CALL  MAL'/E(xAMY,EWHM1) 
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'MrrjN  = M. 

PO  ? f 0 I = 1,N'13P"T 
OTFF~A°S  f FWH  M1  -r  *>H  H ( I , 1)  ) 

Ir  (PiFF.GT.niPMTM)  30  T0  300 

nyrMTNsoicp 
I JK  = ] 

ICC  CONTINUE 

3 F (l  IFNI m.l!  . >0.)  G A TC  301 

rriNT*,«  iNSA:ii"f  cnrv  PEFtREN"-',  fw-m  df  unknown  is 

C"  CHf NNElS  *• 

GO  7 (.  1C 
301  rnNT7*UF 

no  3 l r i=i,nstor 

STuF  EF  (I ) aRFFSTil  c ,t  JK> 

3p?  continue: 

NNX  = c 

7F  ( (NRAi'K.Fo.  0)  .an-).  (*’PKKO.  EO.ni)  30  TD  ’4*0 

IF  (f.^flCK.EQ.a)  I?  70  ’70 

IF  ( r -M 1 > .LT  .01  S’  (1)  ) CRd)  = DT  > T ( 1 ) 

IF  < PP  (N  I)  .GT.  ni 'T  dCHf  NN) ) P9C1J)  = DISTMCHANN) 

3'n  CALL  PC’. in  ' 

NNX  = NNX  * ) 
no  MO  I-! , NCH ANN 
TF(tAMX.I))’»20f*',n,^10 
*,10  MWX(l):  A X Cl ) -nA  C<?  ( I , 

IF<R/*wx<I).lT.0.f>  >AWV(I)  = Qtn 
L2P  CON7  I Mi  IF 
44  (J  CALL  FTr  » FK 

IF  ( t P K :i  .NE  .3)  SO  r3  r 3 0 
4C0  DC  MC  1 = 1, NT 

47C  crxprsn)  = PEK°3S(I)  + DFLTA3 
lc  (LNX)  3'*  0 , 17  0,  3 ?1 
ECO  WFIT  r { ,,  ’o) 

GO  T ( >.*  0 

5 ?0  WR1TF  l >,-0) 

C ...  FAXF  SURE  THAT  Ml  PEAKS  ARE  [ N THE  7 A X G F OF  THE  SPECTRUN 

5 3 3 M = r 

HO  5*  G 1=1, NT 

IF  ( ( »-<FOS(I) -Or_TA?-r.*OcTEN) ..T.OIST(l) > GO  TO  53C 
TF  ( ( •'-KPGS<I>OrirMM.*orTFN>.3r.DIST('OHANN)>  GO  TO  550 
m = v ♦ i 
5FC  CONTI  JLE 

P°  1 M * , " NAOJ  = " , N A 0 I , •*  1 = ”,  1 

IF  (K.(E.NT)  GO  TO  r<qn 
00  5F  0 1 = 1 ,5^ 

560  S T 0 R f r 1 1 ' = J.O 
HA  n J = »/\njn 

IF  (Nt  ",i  ,L  T .10)  GO  TO  L 0 1 

PP I f'T  • UNSAT  ISr  A CT  ORv  PEAKS'* 

GO  TC  *5 
080  NT  = « 

PR1TMM‘.2)  (PFKPOS(T)  ,1  = 1, NT) 

00  f I 1 1=1,7 CHANN 
600  CTOKV.'m  » PACK?(I) 

CALL  ►tSlT* 

K"  PU'.*  =1 
CC  T i «0 
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1C  FOR  '*  T<  - filO) 

1?  FQiflA*  ( 14,412,14  ) 

1*  FO'  «l  » C 1 ^ FG.O) 

If.  FGi\MA  r ( T*  ,13F5. 3> 

lb  roRR/T  ; ) io.3) 

?W  FGi  MAT  (F‘>  .4,11  ,C1|]  .>) 

*0  FORRA  r ( 1H1, 1<  ,M10,10y  , A1C) 

7?  FOrvMfT<10H3  4^8MN8AC'f  = I L , lx  J -H  PKRO  = I2,4X7HNREF 

l,'*yQ>-  ■+X8RNr‘SWl  = I2,*<3HNSS4?  = 12, 4X9HMAXREF  = I 

3u  FORMAT  < 1 1H0  NOHANN  = T',GX7H  ) 

7ri  FOKR/T  CfHO  INIt:AL  G Mr  SSES  OE  pM<  POSITIONS) 

40  FJPMt *<l dHG  CALC  Ji.  A f^O  ESTIMATE"  3F  PEA<  POSITIONS) 

42  FQf?y;  * df  F12.3) 

44  KGrvMA  r (1H0,I4,  It  H °EA» <S  DELETED) 

FNH 


FUNCTION  FACT  O'?  ( “ ) 
nOU?L  E PRECISION  FI, RUM 
11  FACTOR  = 1. 

IF(N-l)tO,*,0,l  3 
IT  IF (N-1P) 21, ?1 , 31 

21  PO  2?  1=2, N 

FI=I 

23  FACTCR  = FCCTCF»  FI 

GO  TO  40 
31  SUM=f. 

DO  3*-  1*11,  M 
FI  = 1 

T4  SUM=EUM*OLOG  (FI) 

T5  FACTCP  = 3L2J3  00  .*OFXP(S'IH) 

40  RETURN 

END 


FUNCTION  °G  A UR  S ( X , A\/  pp  fi  G , SIGN  A ) 
DOUBLE  PRECISION  Z,7SO 
7= (X-AVEPAG ) /SIG-A 
7S0= ( 7**2) / 2 . 

IF(7fO.GE.200.)  GO  TO  10 

PGA U5S=0. 3989422504/SIGMA*OEXP(-7SO> 

GO  TC  11 

1"  PGAURR=0. 

11  PETUPN 
FND 
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SUPPCUTINE  CCN  JJ5(^3cr,H3Ffti<,STn-^,MCHANS|) 

CIMENSION  SPpC  ( *»  C 0 0 ) , n 1ST  (£*100)  , 31 1 5 T (200  ) , I PE AK  (5) 
PIMEtSTON  I HIT  E ( ' ) 

PRINT*,"  NDrAK  = ” » N D E A K , " STT^  = ",STPEV,"  NCHANN  = 
READ*  , (I PEAK (:>!)  , *1=1,  NorAK) 

PR  INT  ♦ ,"  TDEAK  = ”,(TPrAK(N)  , N= 1 , V OF A K) 

PEAP*  , (IHITF (N) , N = l,  1PPAK) 

PRINT*,"  IHITE  = ”,  ( I H IT  E (N ) , N=l,  WZt  K) 

CO  ICC  1=1 , N CHANS 
100  S°EC ( I ) = G • 0 

00  1C00  J=1 , N ° E A 
A VER^  G = IC  EA  K (J ) 

SIGm/'=ST'ifv 
A = 3T.e78  * I HIT  E ( J ) 

00  2C0  I=1,NCHS'4M 
X = I 

SPEC!  R = A*  DG  AUSS  <*,  A '/EP  AG,  SIGMA) 

I F ( I .LT.  AVPRAG  ) GO  TO  101 
IFtSFPCTR.GE.  C.l)  GO  t0  ici 
SPECIF  =0.0 
GO  TO  102 

I'M  IF  (SPFCTR.LT  . 1.0)  SDECTR  = 1.0 
in?  KPECm  = SPEC(I)  * GPPCTR 
M0  CCNT1NUE 
10  0 n CONTINUE 

NSEET  = 200  1 

CALL  PANSET (NSEEH) 

AVEc-rC  = 20. 

IMAX=  2*  Av/  ERAG 
CO  LOO  1=  1 , 1 MA  X 

loo  PIST ( I )=PDQISS (I, A VERAO 
IOLD= 1 
NP=  0 

PO  500  L= 1 , 1 M A X 

A = DIST  (L)  *IMAX 

IF ( A . LT . 1.0)  GO  TO  GOO 

INFW=IOLO*A 

00  Ap0  Jl  = IOLO,I'JrW 

POIST  ( J 1)  = L 

N P = N F + 1 

AGO  CONTI NUF 

IOL  0= I NFW  *1 
GOO  CONTINUE 

CO  6C0  1=1, NCHA NM 
IJ=RPNF(N) *NF 
IF(IJ.LT.l)  GO  TO  600 

SPEC (I )=SPEr (T) *SORT (SPEC (I) )*(30IST ( I J) / A VERAS- 1 . 0 ) 
POO  CONTINUE 
PE TURN 
ENO 
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SUBROUTINE  M AL7E C 3 C , ^WHM^RWHM) 

REAL  LWHM.RWRM 
DIMENSION  S PFC  (210  0) 

P<MAX=100 . 

CO  ICO  1=13 , 2133 

IF  (SPEC  (I)  .LE.OK-iX)  GC  TO  100 

PKMAX=$PEC ( I ) 

IPEA*=I 
130  CONTINUE 

HAFMf X=PKMA X/ 2. 

J=IPFA* 
l-'O  J=J*1 

20  0 IF(SFEC<J) . GT.^A-MAX)  r-0  TO  150 

RWHM=J» (SPEC  (J)-HAF^AX ) / (SPEC( J-l ) -53EC( J) ) 
J=IPFA< 

2*0  J=J-1 

IF(SPEC(J)  .GT.HA- -tax)  C-0  TO  250 

LWHH=  j«.(HArM\X-S~EC(  J ) ) / ( EC  ( J * 1 1 - 5 ° EC  ( J ) ) 

FWHM  = RHH 1 • l_H  PM 

PRIN^,"  = m>'»KHAX,"  LWHM  = " , L WHM , '*  FWHM 

C ,RWH^ 

RETURN 

FNO 

SUBROUTINE  RCS°E:( V, V,M,MAX) 

DIMENSION  Y <4100  ) , X (2F6)  , MW  ( 16  0)  , f (52  0) 

N=  1 
«=1 

10  BUFFER  1N(1,  C>)  (MW(  1)  ,MM160)  ) 

KK=Uh  IT  ( 1) 

CALL  E-RSET  < "COUNT  » 100) 

OECOCE  (150,6  C,MW)  M,NM, y 
IF(KCUUT .GT  .3 ) 30  TO  10 

3 0 FORMAT (2X, I*, 2X , * -A , ’ TFf . 0 ,9 ( /25c 3 . 3 I , / , i c i . 0 ) 
IF(KK.EQ.O)  GO  TO  31 
DO  80  1=1,255 
Y (N)  =X  (I) 

Q0  N=N  ♦ 1 

IF(KK.NF.O)  GO  TO  13 
31  N=N- 1 

ITIMF  = Y(  1) 

IF  ( ITIME  . EO  r C)  rrrMF  ='0000 

PRINT*  ,"  IT  I ME  = , I T T ME 

PRINT  1010, M 

IF (N.LE.MAX)  GO  10  39 

PRINT*,"  N GREATER  THAN  MAX,  N=  ",N 

N=MA  X 

3 9 NN=(N/10>*1 
00  90  <=1,NN 
IN=  ( K- 1)  * 10 
1 1 = I N * 1 
i J=I N*  1 0 

PRINT  1020,  IN,  (Y(  J)  , J = TI,IJ) 

90  CONTINUE 
RETURN 

1010  FORMAT ( 1H0 , T El , "T  A GHORD  NUMBER", 5<  , r 5 , FX, "PFA9") 
1020  FORMAT  ( IX, IF  , 1 CF1Q.0) 

ENO 


SUBROUT  IN E RESIT1 

01  MENS  ION  '»T(M,  W)RH(?q) 

DIMENSION  PEKF'-;R(?01)  , P KBA C K ( 2 1 0 ) , X S EC  T 5 < ? J 1 ) , XSECER(20C), 

1 ORELT  ( 230)  , rX  r,30) 

COMMON  /A  00  / VMM(36?),  IHOL0(2O>,  N 3 3 W 1 , N ]SW2 

COMMON//)  3 1/S  TOR7  f21)‘,,lO),S(210'>>»I7STRT(7C)>I7STOP(23),OJM(3533) 
COMMON  ZftOZ/  SF.-ICG),  ERR7  (20).  *E  ? STO  f 21 ) , 3TOKHH(20)  , 

1 X I NT  IS ( 23 ) , I N 3 r X ( 7 P ) , INITAU23),  IFINA_(20),  STORVl ( 71 0 J)  , 

2 RA1-'X(2103)  , ITTlri,  MNX 

C0MMCN/AC3/  T I Tift  20)  , r EKD0  S (23)  , I r I X pt<  ( 20)  , 3ET A , G ST  A 7 c , 

1 GSOVAL,  NT,  NCC.LC,  ICS,  LO 

COMMGN  /AO-/  INTEEPf  70, 1 0 ) , S 7CRrr < 21 0 0 ) , 0 E .. T A 1 , 3 FL 7 A 2 , D ELT A 3 , 

1 STCIPR,  RFFMAX,  131,  102,  IOT,  NT,  NE,  NRER 
COMMON  / A 0>  / nAC<?(710F),  3R(50',  OACMFO),  NJ,  NOACK,  NPKRO 
COMMON  /OATA  / -3131(2100),  COUNTS  ( 3 1 3 C ) , R J V I ^ , T00AV,  FP.tO,  OFTEN, 
1 NORTON,  NCHANN,  II,  TJ,  TSPECT,  rrMT,  IE RR, MA X 5 EF , MAX° ? 

DATA <F*<T  = 8H ( 1XF9. A ,,  7^3013.2,,  H , 5HE13.2,  3H/4GX,,  oH  9*10.2)) 


OATA 

(woo 

= 

6 m 

3, 

1GX,  s, 

o H 

23X, 

7, 

6H 

30X , 

t,  8 H ♦ 3 X , 5 , 

9M 

5 OX, 

♦ » 

5 OX , 3 

, <JH 

7 0 X, 

2 , 

8 H 

80X, 

1 » 3H  9, 

9 H 

10X, 

9, 

8H 

2 0 X , 7 

, oH 

3 0 X , 

6, 

8H 

40X, 

l> , dM  50 X , 4 , 

9H 

feox, 

T , 

8W 

7 OX,  2 

, PH 

OCX, 

1» 

EH 

9 , 3H  1 3 X , 6 ) 

OATA 

(ISF  = 

IMF),  (MLN<=1H 

) 

• • 

.INITIALISATION.  . . 

JO  = 

0 

JOB  = 

0 

NNX  = 

1 

MTP  = 

1 

MTT  = 

1 

NT  1 = 

NT 

• • 

.STOR 

E INFORMATION  ANO 

GET  RE 

a 3 y 

TO 

CALL  3R 

FPAR. 

• • 

30  00  7C  1=1, NT 
J = NNX  4-1-1 
IHOIP(I)  = IFIX^tJ) 

70  PERSTO(I)  = oR<°73  ( )) 

=0  OElTAl  = 3ELTA1/DRTEN 
OELTf  2 = OELT  A2/3rTEN 
CALL  PRE°AR 
DO  I*-  0 1=1,  NT 
J = NNX  4-1-1 
1.0  IEIXFK(J)  = IHOLDd) 

NNX  = NNX  4-  NT 

IF  (NSSW2.G1  .0)  WUrc  (f  ,90t*2) 

9 0 2 FORMAT ( 1 0 HO  CISTA NCE , 2X8HRAW  QA7A,?X10HFA  C< GROUND ,2X9HG0MP0SITF,  v X 
1 16HRES0LV  EO  SPECTRA) 

n ...  CALCULATE  AREAS  and  EPR0R3  r3R  TMI3  GR3J3... 

00  91  1=1, NT 
XI NT IS  (I)  = 0. 

00  3005  J=1,MAXP? 

30  0 r XINTIS(I)  = X I NT  I 3 ( I ) ♦ ST0R7(J,r» 

J9= JO* 1 

PEKPOS  ( JO ) =P  FR  3 T 0 (T) 

PEKFFR  < JO) =SIGMA ( I) 

PKBACK  (JO)  = 0.0 
XSECT3 (JO) = X INTI  3(1) 

91  XSECER  (JO)  = XSECT3U0)  ERRT(I) 

C ...ADO  3A  CXGROUNOS  AND  PRINT  THE  RESCLVED  ANO  COMPOSITE  SPECTRA 


JPP  = 0 

D03300  I=1,ITIMES 
MOP  = jpo  4 i 
JPP=M0P*IN0EX (I) -1 
M77= INIT4L  (I) 

MMX=1FIN4L ( I) 

IF  <FSSW2.rO.G)  GO  I 0 1530 
IF ( M7Z -M [ T ) 2011, ’Oil , 2010 

2010  MZP  = MZ7  - 1 

tlO  2012  l = MTT,  M’P 

201?  WRITE  {6,rMT  ) (GIST  <L>  , FOUNTS  <L)  ,?10<?  (L) 

2011  WRITF (6 ,93j  2)  <Pr?ST0  ( JP) , JP=M0 3, J P P > 

9C  ? 2 FOFMCT  (40X, 9r 10. ?» 

lr3  0 DO  3 C 10  J-  N 3 ’, MM  X 
K1  = 0 
<2=0 

M = J - M77  4 1 

STORVl(J)  = R4CK2(J> 

DO  1310  JP  = MOP , JPP 
J 2 = JP 

Ic  ( (M.GE.I73TRT < JP) ) . ANO . < M .L E. I ? 5 T 0 P < J3 ) ) ) GO  TO  132J 
1310  <1  = K1  * 1 

IF  (NSSW’.EO.O)  GO  TO  7 0 1 0 

WRITE  <6fpMT  ) OISM  J)  , COUNTS  (J)  ,O10<’(  J) 

GO  TO  3010 

lT’O  DO  1330  JP  - J? , J°P 

IF  ( (M.LT . I7ETRT ( JP) ) .OR. (M. GT .1^ STOP ( JP) ) ) GO  TO  1340 
J9  = JOB  ♦ J0 

PKPACK  ( JP)  = PK3ftG<(JP)  v PACK2U) 

1?’5  <2  = K2  4 l 

ITT  = M - I ’STRT ( JP)  4 l 

STORVl(J)  = ETORVKO  + STORZ  (ITT, J3 ) 

STOR7<ITT,  JP)  = fro?’(ITT,JP)  4 p 1 ? < 2 ( J ) 

1 3T  0 S < <2 ) = STOR?  (ITT, J3) 

1 3U  0 FMT  (3)  = WORD(<l  4 1) 

IF  (ESSWZ.Ef'.O)  GO  TO  ’010 

30  0 9 WRITE  (f,,  FMT  ) ( OI  ST  ( J)  , FOUNTS  ( J)  , 31  0<  2 ( J)  , STORV1  ( J) , (S ( <)  , K = l»  <2) ) 

3010  CONTINUE 

MTT  = MMX  4 1 

3300  CONTINUE 

JOB  = JQ9  4 ST 

02  IF  (NSSW2.-O.0)  GO  TO  ’401 
003400  I=MTT ,NCH1NN 

34G  0 WRITE ( 6 , r MT  ) OI ST (I)  , FOUNTS ( I ) , 31 0 < 2 ( I) 

C ...PRINT  READINGS... 

3411  WRIT  E (b  ,»5 ) (TIT.E(I)  ,T  = 1,9  ),TOOIT 
WRITF (f, 200) 

5-0  DO  5C0  1=1,  JC 
IPF  = IBLN< 

IF  (IFIXP<(1 ).GT.0)  IPr  = 1 3 F 
ARFA  = XSEOTS(I) 4=K3AC<(I) 

STXEPR  = 0. 

IF  ( AREA ,GT .0 . ) >TXERR  = SORTURd) 

WRITE  (fc,  21?)  DEKP0  3 ( I ) , I PF,  P EKE  3R(  I ) , XSF 0 T > < 1 ) , X“ ECFM  I ) , 

1 P K P A C K ( I ) , S T X E ? R 
r90  CONTINUE 

DE VST  = PEVSO  = C.  3 


AUNK  = APES  = 0 . 0 
DO  5 000  IJ<=1,MAXREP 
AUNK=AUNK«-COUNTS(  T JO 
ARES=APES*5T09Vl  (I  JO 

OEVSOsOEVSl ♦ (COUvrS(IJX)-STORVl(IJO  )*»? 

5010  CONTINUE 

DEVSI=SP~T  ( D EV S'' ) 

PRINT*»“  ARES  UNKNOWN  - "jAUNK,"  AREA  CESDLVEO  = ** , A R ES , **  RMS  CR 
COR  = , Dr  VS  T 

RETURN 

=5  FORNAT(1H1,1X,8A10, 10V, A10) 

210  FORMAT (/bX4 HP ESK, ’3XXH  , 19X6HNJM1.  , 3EX+HEST . /A X6MD0SITT0N, S X7H 
1 ,6XGH  ,SX5HPFROR,7XbHCDJjrS,lE<5HEPRDR,13XlGHRAC<GROUNT 

2,6X5HERK0R/ ) 

">r  1 FORMAT  (3X^9 . 3,lXAl,XXPt’.A,4XFa.A,«Xr5.2,5Xr9.2,3(iiXFS.2)) 

2?  2 FORMAT  (Fj.-,lX,?rt0.A,?Fl0.?,F‘?.?,-3.1,3XA:>,<r5.1) 

2?  1 FORMAT (3XF9 .3,1XM,3XFc.4,16XF6.’,3<P9.2,3<4XF3.2>) 

222  FORMAT  <F9.4,lXFlC.4,t0X,2F10.2,Fi.?,-t>.l,TXA6,-3.i) 

?1  ? FORMAT  (3XFJ  . ?,1XA1,‘>6XF6.2,5XF1?.2,3(4XP11.2>> 

■>15  FORMAT  (F9 . ♦ , 1 IX,  2r  1 J . ? , c 8 .2 , F<j  . 1 , 3 X A £ , FE  . 1 ) 

2X0  FORMAT  (2H-t) 

END 


SUBROUTINE  SFCTR* 

C ....READS  THE  S°£CTRJM  CA^OS  A NO  SOS  THE  rrS7ANCr  AND  COUNT  ARRAYS 
DIMENSION  YSAVdD 

COMMON/DATA  /X  (210  J ) , Y (*>10  0)  , RUNT  3,  733  AY,  -RED,  DCL  X,  NOFTEN  , 

1 N,  1 1, 1 J,  I SPEC'S TeMT, IFRR ,MAXR£F,MA<P2 
COMMON  /OPT  N/  <£PR Er , X rPQAT 

C INITIALISE  COUNTS  ARRAY  TO  7ER0 

DO  ICO  1=1, 2100 
100  Y(I)=0.0 
N=  1 0 

DELX  = 1.0 

PRINT*,"  ENrEREn  S°FCTRUM,  CALLUS  RDSPFC" 

CALL  RDSDEC ( Y,N,  NT  AS, MAXREF) 

IF(N.GT.2100>  GO  TO  2 3r 

D1=0.0 

GO  TO  260 

2 *5  WRIT  F (E  , 2 ) M 
25  IERR=1 

GO  TO  300 

D ....  SET  THE  OISTANCE  VALUES  .... 

2r  0 X ( 1 ) = 0 1 

DO  2(0  1=2,  N 

2)0  X(I)=D1  f FLOAT ( I -1)  *DrLX 

3 0 RETURN 

’0  FORMAT (F6.0,9F7.0,6Xp3.0) 

1 FORMAT (6H0  O = ,r3.x,lDH  LESS  THAN  31  = ,r9.3) 

*»  FORMAT  ( 1 HO  NCHANN  = ,T5,13H  GRATER  THAN  2100) 

5 FORMAT (Fo .0/10 (F£ , 0, 2X) ) 

END 
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SUBROUTINE  »EFPr< 

C PETEP  MINES  T3  E PA7A  HE  T "RS  OF  THE  ? rr  ER ENG  F PEAK 

COMMON  /AO’/  OUMftMS  ( tM  ) , ? T ORV1 1 ? U 0 ) , *AWX(213C),  ITIHES,  MNX 
COMMON  /AO 3/  S<I 3 ( ) , L 0 

COMMON  /A0  + / I Nr  rr c(  ?0 , 1 0 ) , S TORr r < 2 1 C 0>,Or:.TAl,nrLTA2,nFLTA3, 

i st oi pr,  refmak,  m,  in?,  iot,  ni,  ne,  n?ef 
COMMON  /A  0 5 / BAG*'?  ( ?10D  , 0n<:>0>,  3ACK(fO),  NJ,  N9ACK,  NPKRP 
COMMON  /OAT  A / OIST(210C>,  C OUNTS ( ? I ] G ) , RJNIH,  TODAY,  FREO,  0F7F 
1 NOFTFN,  N CHANS , II,  TJ,  ISPFCT,  I-MT,  I r RR , MA X * EF , M A X P 2 
COMMON  /O °T N/  KE'RE- , KrPDAT 
IF  (LO)  10,30,10 
10  IF  < NREF)  3 3 0 » 20  C , 300 
*0  IF  (NRFF)  100,’0r,80 
C SELECT  RPF3HAPF  FROM  I NT  rRNAL  DATA 
°0  0090  1=1,20 

TO  STORFF(I)  = INTPFP (I , N'EF) 

C FINO  LOWER  EDGE  Or  3 r A < 

1°  0 00110  1 = 1 , MA  XREF 

IF (STOPEF ( I ) ) 11^, 110,120 

HO  CONTINUE 
1*>0  IMIN  = I 

IF  (I.GT.l)  IMIN  =1-1 
C PINO  UPPER  EDGE  Gc  =>EA< 

I = maxrpf 

1*0  IF  (STOREF ( I ) ) 1*3,140,150 
140  1 = 1-1 

GO  TC  130 
1"0  IMAX  = I 

IF(I.LT.HAXREF)  IMAX  = T«-l 
C FINO  TOP  OF  PEAK 
PEFMf  X = 0. 0 
00  170  t = IM  IN, IMAX 

IF  (F EFMAX-STOREr( I) ) 160,170,1*3 

150  REFMAX  = STOREF  ( I ) 

IT0P=I 

V 0 CONTINUE 

SREF  = .001*  REFMAX 
C FINOS  THIRD  HEIGHT 

REFT  HD  = RE  F MAX/ 3 • 

00  1E0  I = I T OP,I HI  X 
IF  (STOPEF(I)-RE-THO)  190, IPO, 13] 

13  0 CONTINUE 
19  0 RI  = I 

OFFSET  = (RFFTHO-SnREF(I)J/(SnREF(I-l)-STOREF<D) 
RTHPOS=PI-OFFSET 

OFLTA1  = (FL3A'‘(THAX)  --'THPDS)*OrrrS 
OELTA2  = (RTHPUS-rL3AT (IMIN) )*Orr-N 
OELTA3  = (RTHPOS-rL*tAT  (ITOP)  )*OFTEN 
RTHPOS  = (RTHPOS-1 .0) ‘OFTEN 
STORE  = 0.0 

230  I0T=(0ELTA2»CELTA1) 'OF'EN  ♦ 1.1 
101= (OFLTA1 ♦OELTA3) /OF’EN  +0.1 
I02=(0ELTA2-0ELTA3) 'OFTEN  *.l 
NE  = NCHANN  - 191 
N9  * 102  ♦ 1 
IF  (NRFF)  230,210,  230 
f)  FELECT  REFSHAPE  CRDN  RAW  SPECTRUM 


?10  KU  = (STOI^R-OISf (1) ) /nPTEN  ♦ 0.1 
110  = KU  - 102 
IMIN  = 1 

IMAX  = I MI N ♦ IDT  - 1 
STORE  = 0FTEK*FL04r(Iin)  ♦ OISTm 
RTHPCS  = ST  0 iPR  ♦ DELTA  3 
DO  220  IR  = 1, IOT 
??  0 STOPEF(IR)  =RAWX(iin*IR) 

230  IF  (IQ)  24U,2L0,700 

2 0 DO  2*3  I=IMlN,IMK 

2*0  STOREFCI-IMim-1)  = STOrEF(I ) 

DO  21 5 I*=ITOP,IrftX 

2 ' F IF (ST  OREF ( I h > .LT . 3REF ) GTOREF (IM) =3.0 
IM=IMAX-1MIN *2 
00  2t  0 I=IM, kAXRFP 
23 0 STORFF(I) =3. C 

IF  (E  NX)  **30  ,4  00, 3G1 
310  WRIT  E (F  * 310  ) 0ELr4i,  P rLT  A2,  OELTV3 

310  FORMAT  ( 11H1  0 ELT f 1 = , Fy  . 2,  8X9*1=.  r * 2 = , r r . z , a x "mEL  T 4 * = ,F*.2) 
PRINT*,"  IOT  = 'MOT,"  101  = ",131."  102  = ",I02 
WRIT  F (6 , 350  ) RT-I  pO S 

3r  0 FORMAT (32M  REFERS  NOE  PEAK  THIRD  HEIGHT  = ,F9.4> 

410  PETUfN 
END 


FUNCTION  PPOISS( V003, A VERAS) 

1 FPOISS= (( A VERAG* • NO 3S ) /F ACT  OR < NOT S > ) * EXF ( - A VERAG > 

PETIJCN 
END 


SURPCUTINE  OUTP'lT(X,NMAX) 

OIMR  SION  X ( 1) 

IMAX=NMAX/1C 

00  ICO  K= 1 » IM  AX 

1 1 =K*  10 

10=11-9 

inn  FRINT  200, (X (I) , 1= 10,  Tl) 

?10  FORMAT (EX, 10  (E9.3,  2X)  ) 
PETURN 
END 
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SUBROUTINE  A II 


SUBROUTINE  PjRNn 

C , • .CALCUL AT  ES  THF  1 \ C Kr-ROUNTS  • • • 

COMMON  /A  02/  Df(;?),  P“02(:2>,  0 5 0 * * X < 1 0 > * 4 ( 1 C , 1 0 ) , J'JMMY  (.  04  0 > 

COMMON  /AOT/  TIT  „F  (20  ) ,PEAK(  20)  , 3.X<  ( 23)  , X 3 r A < , NT AlC , I3S , 10 
COMMON  /4  0* / INTarrf?0,10),STOREr<?100),3F»Tftl,OELTA?,JELT4?, 
i ST01PR,  REFMAX,  111,  102*  InT,  XT,  HE,  X?EF 
COMMON  /AO,/  1AC<2 (210*1  , 2B(60),  TACK(T3),  NJ,  NBmCK,  NPKRO 
COMMON  / HA  T A / r»1 5T  ( ’ 1 0 0 ) , COUNT*  < ’13  0 ) , RJXID,  TOOAY,  c3.cr' , 3^7ci, 
1 NOFTFN,  NCH4NN,  S<  I ° < 5 1 , **  A XREr,  14  < 3 2 
00  = OIST(l) 

IF  (NBACK-l)  10,TQn,9GP 
r ...AUTOMATIC  T4CKGR0UN0  SELFCTnX  . . . 

10  NP  = NPEAK  ♦ 1 
IJ  * 0 
NJ  = 1 

RH02  (NJ)  = 1 .E10 
00  2C0  JJ=1,NP 
IJ  * IJ  M 
IF(IJ-NP)  20  ,3  0,-  0 0 

?0  LOOPY  = ( ?F  A K(  I J)  - 0ELTA3-Q”')  /OFT  EX  ♦ 1.1 

IF(IJ-l)  40,40,60 
♦ 0 NOEL  = LOC3K  - 102  - 1 
MA  = 0 
GO  TO  103 

' 0 NOEL  = 1FIX  ( (PEAMT  J) -PEAK(I  J-1D  OrT  EN  ♦ 0.1)  - IOT 
MA  = LCCP<  - 102  - XOPL  - 1 
GO  TO  IOC 

10  LOOPY  = (P£AX(IJ-1 ) -1FLTA3-QO) /O-TF X ♦ 1.1 
NOEL  = NCHA  4N  - _0C3K  - ID1 
MA  = LOC  = < ♦ ID1 
100  IF(NCEL-IO)  200, 123, 12F 
120  J8ACY  = FLOAT (NOrL)/10.0  ♦ 0.2 

jo  = ndel/jback 

SHIFT  = FLOAT (JO) / 2 , 

no  lfcO  I=l*JpACK 
NJ  = NJ  ♦ 1 
B =•  0.0 

00  HO  10  = 1, JP 
1*  0 B = ° ♦ CO'JN  IS  ( Ml  ♦ I O) 

RH02 ( N J)  = o/FLOtT ( Jp) 

BACK (NJ)  = RHO? ( X J ) 

CHECK  = RH02  (NJ-1 ) ♦ 2.*SORT(RMO?(XJ-l)+l.) 

IF  (RH02(NJ)  -CHECK)  1F.J,  160,150 
IT  0 NJ  = NJ  - 1 
GO  TO  170 

IsO  OE(NJ)  = OIST(MAH)  ♦ OFT  EN*SHICT 
09 (NJ)  = ( NJ) 

C ...TEST  FOR  LIMIT  Or  b 0 POINT3  AXO  Tt  <E  APPROPRIATE  ACTION... 

IT  0 IF  (NJ-,9)  180,1*3, ,00 

lfF  If  (IJ.GF.NP)  GO  TO  A 0 0 
IJ  = NP 
GO  TO  80 

13  0 MA  = MA  ♦ JP 
200  CONTINUE 

ADO  IF  (NJ-6)  *10,420,420 
♦10  M * NJ  - 1 
GO  TO  430 


V 
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u?n  m = l 

'*■*0  N = NJ 
500  MM  s H ♦ 1 
00  510  1 = 1 1 M 
*510  Mil)  = 1.0 

CALL  LSQPQUN.MM) 

DO  5L  C 1=1 , NCHANN 
RACK?  (I)  = X(l) 

DO  5?0  1P=1,1 

5?0  BACK  ? ( I ) = BACK?  ( I ) ♦ Y ( I P*  1 ) * DT  3 T ( I>  **I  0 
IF(BACK?(I> ) 530,5  43,5*  0 
R30  BACK ? ( I ) = 0.0 
0 CONTINUE 
GO  TO  800 

C ..  .INTERPOLATE  HE  °ACKGR  OUN")  THROUGH  THE  PJTNTS  REAP  IN... 

■’0  0 JC1=C 
JC2=  0 

IF(OB(l)-OIST(t) ) 710 ,7 10,?05 
■’05  JC1=1 

OE ( 1 ) = 0IST  (1) 

RHO? ( 1 ) =BAC  K (1) 

710  Ir(0R(NJ) -HIST (NCHANN ) ) 715,750,753 
71E  JC  2=  1 

OE(NJ»JClU)  = 0 T 3 T ( NCHANN) 

RHO?  ( N J*  1C1+  1)  =Bf.  CK(NJ) 

’■0  no  715  1=1,  NJ 
OE  (X*  JC1)  = BB ( I) 

7-G  PH02  < I J Cl ) = RAC<  ( I ) 

N=NJ ♦ JCi*  JO  ? 

IF  (N.LT.3)  GO  TO  910 
J=  2 

00  790  1=1, NCHANN 
OIX  = 0I5TCI) 

7»0  IF ( DI X -OE ( J M > ) '50,77-, 775 
7*5  J*J»1 

IF  (J.GE.N)  J=  J - 1 

7 8 0 OEMl=(OE< J-1)-3E< J)  > * (nE ( J-l ) -OF( J + l ) ) 

DEM?  = ( OE ( J ) -CE(  J-i  >>  MOE(  J)  -DP(  JM>  > 

0EM3= ( D E(  J M > -BE (J-l) )’  ( OE ( J *1 ) -1 E ( J ) ) 

BACK?<l)=('HX-OE(J>>*  <nIX-PE  ( JM)  1 / DEMI*  RH  02  (J-l) 

1 ♦ (OIX-OE (J-l) ) * (niX-OE( J*l) 1 F3EH?*RH3?( J) 

2 ♦ (OIX-OE  (J-l ) ) ‘ ( 01 v-DE( J ) ) /OEHT*  S HO? ( J ♦ 1 ) 

IF(BACK?(I) .LT.0.0)  BACK2(I)=0.1 

7 ■*  0 CONTINUE 

800  WRITF(6,313 ) (Or(I),RHP2(I) , 1=1. N* 

810  FORMAT  (?rH3POINT“  JB^O  BY  BA  CKGO  JM5  / ( ?F  1 2, 3)  > 

GO  TO  1000 

c ...calcjlate  the  background  intervale... 

BOO  IJ  * 2 

OB  ( 1 ) = OIST(l) 

IF  (NJ.EQ.l)  OB  ( 2)  = 0TST  ( NCHANN)  0 FTF  N 
DO  9?  0 I=1»NCH4NN 

IF  (DIST  (I)  .GE.OBdJ)  ) IJ  = IJ  * 1 
IF  (IJ.EQ.NJ)  OB(MJ*l>  = OISKNCHSHH)  ♦ JrTEN 
9?0  BACK?  ( I)  = BACKdJ-l) 

WRITE (6,930) 

930  FORMAT  ( 13H0  0 IS  T ANCr  , 17X 1 1H  B\  ?.<i  ROUND) 
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00  9‘  0 1 = 1,  NJ 
DIX  - 0P(I*1  > - :?TrW 
9'.0  W9ITF(6,l*53)  (DB(r),DIY,  BACK  ( I > ) 
9' 0 FORMAT (F9. 3, 1H-,F9 .3, 1* XF6.2) 
1000  RETUFN 
FNO 


SUBROUTINE  L SQPOw  ( NS'.in  , MSU*) 

LEAST  SQ'JA-<F  POLYNOMIAL  FIT 

COMMON  /AO?/  X(5?>,  Vlr?),  W(50),  3(10),  A(1G,10),  XPOWER(tiO), 
1 OUMMY<40h9> 

NsNSOB 

M*MSUP 

M1=M*1 

M3  = M + M +M 

M31=M3-1 

M41*M31*N 

FORMATION  AND  INVERSION  OF  SYST^H  0-  NORMAL  EQUATIONS 

00  100  K2  = M1  ,M4i 
XPOWFR («2) = 0 .0 
100  CONTINUE 

DO  200  K1=1,N 
TERMsW(Kl) 

00  2 C 0 K2=M1,M31 

XPOWFR (K2)=T ERM+XPOWER ( <2) 

TERM=X  (Kl>  * TERM 
200  CONTINUE 

00  3C0  1*1, M 
00  300  J=1,M 
K2  = I ♦ J+M-t 
A ( I , J)  =XPOWE  MK2) 

300  CONTINUE 

DO  4 C 0 K=1,N 
TERM  = W (K)  *Y  ( K) 

00  400  K2=M3,M41 

XPOWFR (K2> =T  ERM** POWER (K2) 

TERM  = X (K)  *7  ERM 
41 0 CONTINUE 

00  ECO  1*1,  M 

K?=I*M31 

9(1) =XPOWER< K2) 

FOO  CONTINUE 

CALL  MATINlM  A,  M,3,  1 , 9ETERM,  1 0) 

730  CONTINUE 
RETURN 
END 
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SUBRLUTI NE  PREFA? 

C ...CONVERTS  CATA  TO  CHANNEL  NUHA'?3  AND  3 RfPA  Rr$  IT  FOR  THE 

C VARIABLE-METRIC  PACKAGE... 

PEAL  LPF 

COMMCN  /A  00/  HH(?0 , 20) ,PER<  20) ,GR( 20)  ,S( 20)  ,XP( 2 0) ,GP (20) ,T (20) , 

1 GB(23),eD,GS,c„,S..,Fo,GSP*T0*',7?,3,AA,SSS,P0»GrP,FD,GTT,GSP, 

2 DELTA,EE,  LT  , MS  , I T , L » 7 HOLE  ( 20  ) , N53H1  , NSSH? 

COMMON/A  01/7  <213  7,n),''S<2102),'',’5TRT<2C),IZSnP<2r!),R<213  2), 

1 X ( 21 00)  ,'4(2130)  , E ( 2 1 1 0 ) , C ( 2 0 * 20)  , A(2n,  GCu),  ERR  ( 20)  , Dp  L > 

2 REP,  F,  P B FST , _P-<22),  LP(22),  I <5 , I3R,  I°RE,  IlC,  ITC,  IOFF 
COMMON  /SO?/  SIO-A(?0),  FRPZ(20>,  PERSIC (21),  ST  ERhh (20)  , 

1 XlNTIS(20)  , IN 0 E X ( 2 0 ) , INITAK20),  IF1NAl<20>,  SI ORV1 ( 21 0 0) , 

2 RA WX ( 213G ) , ITMp3,  MNX 

COMMON  /A 0 3 / DUMMY ( 3 T ) ,NT ,NCALC» 135, LO 

COMMON  /AO,/  INlPrF(?0,10),STORFr(2lC0),TELTAl,nELTA2,DELTA3, 

1 ST  Cl  PR*  RFFMAX,  111,  132,  IDT,  >13,  Np,  NRSF 
COMMCN  /DATA/  DI5T{?10C),  COUNTSt 21 3 0 ) , RJNIB,  TODAY,  PRPQ,  OFTEN, 
1 SK1P(7)  ,MAXRPr,MAXP? 

PEL  = 0.1 
QQ  = DIST(l) 

C ...CONCERT  AL.  UNITE  TO  CHANNE.  NJMPp*... 

00  1C60  1=1, NT 

PERSTO(I)  =(  PERSTK  T)  -On) /OFTEN  * L.O 
tOeO  STCRHH(I)  = 1.0 

PERSTO  (NT*1)=3C0P.P 

DEL 3 =AINT  (DFLT?  1>DELTA2»0. 1) 

IPF=  (2.0OELTA2) 

IPR£  = ( <2.3*DPLTA?)  --  L OAT ( I®R ) ) * 1 0 . 0 
DO  lCbl  I = 1 , M A X 3 2 
10*1  R ( I ) = 0 . 0 

00  1 T65  1=1, IDT 
10  E E R (1*1 ) = STORpF(T) 

IRS  = IDT  » 1 
I1C=1 
ITIMFS=0 
III  = 1 

C III  = (NUMBER  Oc  PEAKS  ALREADY  ANALYZED)  ♦ t 

2000  IOFF  = III  - 1 

C ...DIVIDE  GROUP  Op  CEA  KS  INTO  NUN-OVERLAPPING  SUBGROUPS... 

MS  = 0 
LT  = 0 

IF(III-NT)  ?0  Cl,  ?"01,  3000 

C I1CSS  = (FIRST  CHANNEL  NUMBER  DF  FIRST  PEAK  IN  SUBGROUP ) - 1 

?0n  I1CSS  = PERSTO(ITI)  - PpLT  A 2 ♦ 3.1 
IFdlCSS.GT  . l)  I1DS3  = I1CSS  - 1 
DO  2C 15  JJ= 1,20 
LT*LT*1 

IF  (PFKSTD  ( III* /J) -PF^STO  (IIIOJ-1)-  DELS)  2 C 15 , 20 15  , ? J20 
’015  CONTINUE 
20,0  XXX  = PERSTO(III) 

C ITC  * TOTAL  NJMBER  PF  CHANNELS  IN  THE  SUBGROUP 

ITC  = PERST  0 (LT ♦TDPP)  - XXX  ♦ OE.TAl  ♦ TELTA2  ♦ 3.1 
C ...SETUP  INITIAL  MATRIX,  DATA,  4EIGH1S,  AND  PARAMETERS... 

DO  2025  1=1, LT 
PO  2025  J= 1 , LT 
2025  HH  ( I , J)  = 0. 0 
OELTAM.3 
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00  2 C 3 0 1 = 1, IT 
IM  = I0PF  ♦ I 

PFK  ( 7 ) = PfcRSTO(IN)  - XXX  ♦ 0EL“  1 2 2. 

MH(I,I)=Sn?FH  ( M)  **  2 
IF  ( IH0L0  (I’D  ) 20  26 , 2G  **8 , 20  27 
TO’!  HH(I,I)  = a.u 
GO  Tr  2030 

?0“0  DELTA  = DEL  T A *HH ( I , T ) 

21TG  CONTINUE 

no  2 c 35  i=iic,it: 

JJ  = I ♦ IlCSS  - 1 
X(I) = RAWX ( )J) 

207E  H<I)  = 1.0/SORT(*<Ii  *10.0) 

EE  = FLOAT ( I TC -L  T) 0r 
CALL  DAVION 
ITIMFS  = ITI>1FS»1 

INPEX  = N LM9E  R Oc  P"AKS  IN  SUBSRDlte 
IK'ITAL  = STARTING  CuANNrL  N'JRSrR  IN  E J3SROU3 
IFINAL  = FINAL  CHANNEL  NUM9Ef'  rO  R SUPSRDUP 
INDEX ( ITIMES? =LT 
INIT  AL ( IT  I N E S)  = UTSS 
IFINAL  (ITHFS)  = tlDSS  ♦ ITC  - 1 
03  = ITC  - LT 

C ...CALC.  ERRORS  IN  »nS7TI0NS  A N1"*  ODiPUTF  NORMALISED  INJIV.  °EAKS 

00  2C4  2 I--l,l.T 
IN  = IOFF  ♦ i 

PEFSTO  ( IN ) = (°E? ( I)  *X*X-0ELTA?-3.Gi  * OFT E N «•  00 

STORHH(lM) =HH(I, I) 

SIGMA  ( IN)  =S0  P T (*-PS  ( ( 2. 0*5T  ORHM  ( 1 1 ) * FP"3  n /QT)  ) * OFT  £N 
IF  ( A ( I ) ) 211E*2110,211r 
2110  ERh  7 ( I M)  = 0.0 
GO  TC  2120 

21 1 S FRR7  (IM)  =ERP  (I)  / A ( I ) 

Rl^O  00  2 04  0 J=l,  f*A X° 2 
?n;fl  7(J,IM)  = 7 ( J,IM) • A (I ) 

20 2 CONTINUF 

IF  (IT.GE.2E)  MRITE(6,  1001)  (p  ERS  TO  ( I ) , I = 1 , i.T ) 

III  = III  ♦ LT 
GO  TC  2000 
WO  RETUFN 

1001  FORMAT  ( 51H  ] EPSILON  TOO  SMALL  - OlNOT  CONVERGE  FOR  PEAKS  AT  ,:F10 
1 . 2/ ( 1 0F10 . 2) ) 

END 
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SUBROUTINE:  FCN(L.TT,3F,F9,pER,M1» 

C ...SETS  UP  THE  EQUATIONS  AMD  DETAINS  RELATIVE  PEAK  HEIGHTS 

REAL  LP* 

DIMENSION  PFR(20) , G?  ( 2*) 

COMMCN/AOl/ ?( 210?, i)>  ,SS(210?)  ,I7STRT  (20  , I7STQP  (20)  ,R(2102)  , 

1 X(21Q0)»W(?t00) ,r(?l?0) ,0(20,20) , A ( 2G  ) , G ( 20 ) , £RR<20),  TEL , 

2 REF,  F,  r BEST , .»-(22),  L®(22>,  IRS,  I3*,  IPRr,  IlC,  ITC,  IOFF 
COMMON  /DATA/  QIST(?100>,  COUNTS  (2100)  , RJNTD,  TODAY,  FREQ,  OFTE  J 

1 NOFTEN,  NO-ANN,  II,  7J,  ISPeCT,  I-MT,  I ERR , Ml  X REF  , **  AXP  2 
LT=LL  TT 

IF (Ml . EQ. t)  F9ES  T = 1.0E10 
4006  DO  91  L = 1 , L T 

LP(D  = FER(L) 

LPF(L)  = (PE R( L)  -r  LQ A T (LP(L>  ) ) * 10 • 3 
PI  CALL  SHIFT ( L) 

PRINT*,"  IF  M = 3 : ".Ml,"  WILL  C1_L  FUN  FROM  FON  " 

205  IF  ( Ml  -3)  222,  21?, 22? 

C OUMMY  CALL  TD  FUN  TO  FORM  THE  ARRAY  E 

219  CALL  FUN(LT) 

DO  220  1= II C ,ITC 

270  W(I)  = 1 . 0 / S CRT  (&3S(E(I))+10.) 

272  DO  210  M=1,lT 
MM  = M ♦ IQFF 
DO  210  L = 1 , M 
LL  = L ♦ IOFF 
C(M,L)  = 0.0 
DO  200  I = I1C  ,ITC 

IF  ( (I.LT.I?STRT(LL) ) .OR. (I.GT.I7STDP (LL) ) ) GO  TO  200 
IF  ( (I.LT.IZST.RKMMH  .OR.  (I.GT.r^5TJP  (MM))  > GO  TO  230 
IL  = I - I7  S 7RT  ( L L ) ♦ 1 
IM  = I -IZSTFT(MM)  ♦ 1 

C(M,L)  = C(M,L)  «•  H(I)‘'H(I)»Z(IL,LL)  *Z(IM,MM> 

2^0  CONTINUE 

21  0 C < L , M ) = C ( M ,L  ) 

00  201  M=1 , LT 
MM  = M ♦ IOFF 
A (M)  = 0.0 
DO  201  1=1 1 C , ITC 

IF  ( (I.LT.I TSTRT(MM)  ) .OR. (I .GT.I7>TQP (MM) ) ) GO  TO  201 
IM  = I - I7STRT( AM)  ♦ 1 
A ( M)  = A ( M ) ♦ W(I)  M(I)*X(I)*7(IM,M1) 

201  CONTINUE 

CALL  MATINV(C,LT ,1 ,1, DETER*.  5) 

CALL  DERI VP ( LT  ,°FR) 

00  If  0 0 L=1  ,L  T 
1600  GR(L>  = G(L ) 

PRINT*,"  CALLING  fij^  tpqm  fCN,  c - ",F,”  GR  = 3=  ",G 
CALL  FUN(LT) 

I F ( M J - 3 ) 206,1011,206 
1011  F7  = F/(ITC-IlC-.n 
DO  1C 06  L = 1,LT 
REP  = F7*C(L  ,L ) 

1006  ERR(L)  = SQRT( A9?(REP) ) 

206  F9  = F 

IF  (F9.LT.FBEST)  e BEST  = F9 

RETUPN 

END 
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SUBROUTINE  OERI V = ( LL T T , PER) 

• • .CALCUL A TES  OF^IVITIVES  OF  FUNCTION  F(FJN)  WITH  CHANGE 


IN  POSITION 
REAL  LPF, LOPE 
DIMELSION  PFR(20> 

C0MHCN/AC1/7  (2102,  11)  ,CS<210  2)  ,I?STRT  (20  , I Z STOP (20) , R ' 7 1 0 2) , 

1 X(21J0)  ,W  (2110)  , E( '*100)  ,C  ( 20,20)  , A ( 2 0 ) , G(2d),  ERR(20>,  OEL, 

2 REP,  F,  c BEST , _p-(2?>,  ° ( 2 2 > , I*  S,  I3*,  IPRE,  IlC,  ITC,  IOFf 

COMMON  /DATA/  1151(2100),  COUNTS ( 210 0 ) , R JN ID » T DO AY , FRFD,  OFTEN 

1 NOFTFN,  MCt-ANN,  II,  IJ,  ISPECT,  I?MT,  IERR, MAX ? EF, MA XP2 
LT=L ITT 
00  9fL=l,LT 
LL  = L ♦ IDFc 
DO  9<-  1 = 1,  MA  XP? 

3L  SS(I)  = 7(1, LL) 

ISA  VI  = IZSTPT (L-) 

ISA V?  = I ZSTD  P ( L . ) 

LOP=LP(L) 

LOPE  = L PF  (L) 

Y = PEF-  (L)  + 9EL 
LP  ( 2 1 ) = Y 

LPF ( 21 ) = ( Y-FLDAT (LP (71) )) *10.0 
Y=PER (L) -DEL 
LP ( 2? ) = Y 

LPF ( 22 ) = ( Y-FLOAT (LP (22) ) ) * 10.0 
LP(D=LP(21) 

LPF(L)=LPF( 21) 

CALL  SHIFT ll  ) 

PRINT*  ,**  CALLING  ^UN  FPOM  0ERIV°,  L = “,l,“  F = ",F 
CALL  FUN(LT) 

FP  = F 

LP (L) =LP( 2? ) 

LPF(L)=LPF(  2?) 

CALL  SHIFT (l  ) 

PRINT*, **  CALLING  FPOM  DERIVO,  ra  = F = **,F 

CALL  FUN(LT) 

G ( L)  = ( FP-F)  /<  2.0'DEL) 

00  9E  1 = 1 , M A XP  2 
95  7(1, LL)  = SS  (I) 

IZSTFT(LL)  = ISA  VI 
IZSTOP(LL)  = ISAV2 
LP (L) =LO° 

?6  LPF(L)=LOPE 
RETURN 
END 
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SUBF  OUTINE  SHim.O) 

...CONTROLS  INTE^OLAT^  ON  OF  FRF.  a = AK  FDR  SHIFTING  FRACTIONS 
OF  CHANNELS 

...CONTROLS  IN  TEr  30L AT  TO  N OF  REr.  PEAK  FDR  SHIFTING  FRACTIONS  OF  C 
REAL  LFF 

COMMON/A 3 1/7  ( 21 1 2,  11 ) , CS ( 21 0 2) , I Z 5 TR T ( 20 ) , I ?STD° ( 2 0 ) , R ( 2 1 3 2) , 

1 X ( ? 1 C 0 ) , W ( Zl 00  ) , E(,lnO)  , 0 ( 20 , TO ) , 4(00,  G(23),  FRF(2Q>,  OEL, 

2 REF,  F,  -PEST,  _P-(2?),  L?(22>,  IRS,  IPR,  IPRE,  I1C,  ITC,  IOFF 
COMMON  /DATA/  DIST(TlOO),  COUNTS (? 13 0 ) , RJNID,  'OOAY,  FRiO,  OFT E 1 , 

1 NOFTFN,  NCHANN,  II,  IJ,  ISPFCT,  I-MT,  IErP, M4XREF,MAXP2 
L = LO 

LL  * L ♦ IOFF 
C1  = LF (L) 

C2=LPF (L) 

U=C1*C2/13. 0 
P8=IFR 
°9=IFRF 

00  ICO  J=J.  ,MAXF2 
ICO  7 ( J, LL ) = J.O 

P1  = BP  *09/13 .0 
IF  (FI  -U)  TOO,  IS, 3 
r,  SUBROUTINE  SHIFTL  L LE^S  THAN  R 

9 J=1 

I7STFT  (LL ) = 1 
IF  (IPRE-L°F(L) ) 10,11,12 

10  F = 1 0 *IDRE-L  PF( t ) 

1 = IPR-LP(L) 

GO  TC  14 

11  F=O.C 

GO  TO  13 

12  F = IPRE-LPF(L) 

131=  IPR-LR(L) *1 

14  FI  = F/10.0 

15  7 ( J,  LL  ) = R ( I)  ♦ Fi*(o  ( I ♦ l)  -R(I)  ) 

1 = 1*1 

J=J*1 

IF  ( 1 -IRS)  IS, 15, 17 

16  IF  { J-ITC ) 15, IF, 17 

17  I7ST CP (LL  > = J - 1 
GO  TC  23 

C SUBROUTINE  S*IFTD  L FOHAL  r 

18  DO  20  J=1 , I RS 
20  7 ( J , LL ) = R ( J ) 

IZSTFT(LL)  = 1 
IZSTCP(LL)  = IRS 
GO  TO  23 

C SUBROUTINE  SHIFTK  R LE^s  THAN  L 

3?  0 1 = 1 

IF  (IPRE-LPF  (L)  ) 1,2,? 

1 F = 10*I3RE-LPF(u) 

K = L P (L  > -I  PR*  2 
GO  TO  5 

2 F=  0 • 0 
GO  TO  4 

3 F = IPRE-LPF  (L) 

4 K =LP  (L)-I°R*1 

5 FI  s F/10.0 
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IZSTRT(LL)  = < 

J = 1 

6 7 ( J,  L L ) = R(I)  ♦ :1MRI1U)-R(l)l 
1*1*1 
J=J*1 
K = K ♦ 1 
IF  ( I -IRS ) 7 ,7,3 
T IF  (F-ITO)  b ,6 , 8 
8 I 7ST0P (LL ) = < - 1 

?■*  return 

ENO 


SUBROUTINE  FUN(LD 

...CALCULATES  T-tr  FUNCTION  F WHIH  IS  BEING 
REAL  LPF 

COHMCN/AOl/7  ( 713  2,  1 3 ) , '‘S  ( 21 0 2)  , I 7 5 T R T ( 2D  , I^STnP  (20)  , R( 2102) , 

1 X(  2100  , W (21J0)  , " (21,,0)  ,2  ( 20  ,’C)  , A(2P),  G(2J),  ERM2Q),  DEL, 

2 REP,  F,  F BEST,  *3r(20,  LP(22),  IRS,  PR,  IPRE,  110,  ITC,  I0C 
COMMON  /DATA/  DI“T(?10n),  COUNT S ( 2 1 0 0 ) , R JNI D,  TOOmY,  FREQ,  OFT 

1 NOFTEN,  NCHANN,  II,  TJ,  ISPECT,  I-MT,  IERR , MA X REF  , M A X= 2 
F=  0 . C 

00  60  1=1, ITC 
'0  E(I)  = 0.0 
DO  70  L = 1 , L T 

11  = IZSTRT  ( LdOr:r  ) 

12  = 1 7ST03  (L  *1  OFF) 

C-  0 1=11,12 
J=I-Il+l 

20  E(I)  = Ed)  ♦ A(_) *7< J,L*IOFF) 

00  80  1= I1C , ITC 

FI  = (X(I)-E(I)) 1 W(I) 

30  F = F ♦ FI*  FI 
RETURN 
ENO 


SUBROUTINE  A RYOO R( A,B,N,L) 

ORDERS  THE  AR?AY  A,  AND  MAKES  SHE  ORDERING  TO  S IF  L = 2. 
INTEGER  B 

DIMENSION  A ( 1)  , B ( 1 ) 

NP  = N - 1 
00  5 C 1 = 1, NP 
K = I ♦ 1 
DO  5 f J = <, N 
IF  (A(I)-A(J))  50,50,20 
20  IF  (L-l)  40,^0,30 
TO  IS A V = B ( I ) 

R (I ) = B ( J) 

B(J)  = ISAV 
A 0 SAVE  = A ( I ) 

A(I)  = A ( J ) 

A ( J)  = SAVE 
"0  CONTINUE 
RETURN 
END 
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SUBROUTINE  M AT  IN  V ( ft  , NS'  IB  , 9 ,M SUR , 0 i T , NM  AX  ) 

dimension  a (nmax.ns  n>  ,B<NMAx,M5.m 

DIMENSION  PI V0T( 13  0)  , INDEX ( 1 0 0 ) . I 3 [ V OT  ( 1 0 0 ) 

EQUIVALENCE  ( p IV  D T , I s D "X  , IP  I VOT ) , ( A 1 1 X , T ) , ( IRON,  il.IRO  , 
1 (T EM  = , SNAO) 

OATA  (ISHlFT=409i)  , (MASK=OOJ00003777777775) 

N=NSt'P 

M=MSU8 

INITIALI7ATICN 

0£TFFM=1. 0 
00  2 C 1 = 1,') 

IPIVCT  (I)  = 0 
*»0  CONTINUE 

00  SF  0 1=1, N 

SEARCH  FOR  PIVOT  ELEMENT 

AHAX=0. 

DO  1C5  J= 1 , N 

TEMP=IPIVOT  (J)  .ASH, ,k)07  • MASK 
IF  (TEMP)  1 J 5>  ,60. 

FO  DO  ICO  <= 1 , N 

TEMP  = IFIVOT  ( X)  .A  N^  . . M0T  .MASK 
IF(TEMF)  100,30 
* 0 TEMP=A9S  (A  ( J,<)  » 

IF (TEMF-AMA  X ) IOC, 3 5 , 8 f 
IE  IROW=J 
ICOLUM=< 

AHAX=TEMP 
100  CONTINUE 
105  CONTINUE 

INDEX  (I) =IN0  EX(I)  (ISHIFT*IROH+i:Ou  J*  ) 

J=IRCW 

AMAX  = A ( J,  ICOLUH) 

OETEFM  = AMAX*  CETERA 

MATRIX  SINGULAR 

IF  ( 0 FT  ERM)  110,500 

110  PIVOT  (ICOL'JM)=INOEX(ICPLUM>  . OR . A HA  X . A NO. . NOT . MA S * 

INTERCHANGE  ROWS 

IF(IFOW-ICQLJM)  1+0,  ?F  0 
t'+0  DETEFM  = -nETERH 
DO  200  K=  1 , N 
SWAP  = A (J,  <) 

A(J,K)=A<ICOLUM,<> 

A ( ICC  L UM , <) =SWAP 
pro  CONTINUE 

DO  210  K=1,M 
SWAP=B(J,<) 
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DEVST=PEVS3=  C.  3 
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— ...  


9 ( J,  *<)  =9(  ICOLU*,  <) 

9 ( ICGL  UM , K)  = SWA° 

27  0 CONTINUE 

01  VICE  PIVOT  ROW 

23  0 K=ICCLUM 

AdCOLUM,  <)  = 1.3 
00  3l  0 K=l,  N 

A <ICCLUM,<>  = A<  IC0LU-1,*O  /AMAX 
77  0 CONTINUE 

00  37  0 K=1,M 

9 ( ICCLUM,  <)  =9  UCOLU1!,  K)  / AMAX 
3’  0 CONTINUE 

REDUCE 

DO  5E0  J=1,N 
IF  <J-IC0L'JM)  -» 0 0 , .5  0 
430  T = A ( J,  ICOL'JM) 

A ( Jj ICOLUW) =0.0 
00  4 f 0 K=1,N 

A(J,K')=A(J,K)-A(I3  0LJM,K)+T 
430  CONTINUE 

DO  500  K = 1 1 M 

8(J,K)=B(J,  K)-T(i:0L'lM,<)*T 
500  CONTINUE 
«?r0  CONTINUE 

INTERCHANGE 

500  DO  710  1=1, N 
I1=N+1-I 

IRC= INDEX (II ) .AND. HARK 
K=IRC /ISHI^T 
ICOLUM=I»:-K*ISHrrT 
IF  ( K - 1 C 0L  'J  M ) 650,710 
f>-0  DO  7 CD  J=  1 » N 
SWAP= A ( J, <) 

A ( J,  *0  =A(  J,  1 COLD *•) 

A ( J, ICOLUM) = SWA0 

rn?  continue 

710  CONTINUE 

7'*  0 OET=OETERM 
C RETURN 

ENTRY  DMATINV 

RETUFN 

END 
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SU8RC  UTINE  DAVIS') 

C ...THIS  IS  THE  CONTROL  ROUTINE  *3  ? THE  V A RI  A RLr- k‘ FTRT  D PACKAGE. 

COMMON  /ft  0 3 / H ( 2 r , 21 ) , v (20)  , G(20)  , SC  2 0 ) , <=(  20) , G*(2C)  , T<  23)  , 

1 G8(20),F,GS»EL»SLfcP»GSP»Tn»Z.T»ft»SSS>FO,GTP.FF,GTT,'>^9. 

2 OELTA,E,N, MS  * I T , L,IHOLC(20) , N*SW1 , NSSW2 

1*  Ml  = 1 

IT  = 0 
F=0.  C 

CALL  FCN  ( N » G ,F,X,M1) 

IF(NFSWl) 20 , 25, ?C 

■»0  WRITE  (6,3)  IT, MS,* 

GO  TO  6A0 
120  Ml=2 
200  CALL  READY 

PRINT*,"  CAL  - rQ  R*ADY  L=  -,L 

GO  TO  (800, 3 CO, 500  , *0  0) , L 
300  CALL  AIM 

PRINT*,"  CAL-EO  AIM  L=  ",L 

GO  TO  ( A0 0,f 00,500, >00),  L 
400  CALL  fir* 

GO  TC  (500,500, 3T0, 3 0 P > , L 
*00  CALL  ORES3 

IF  < I T -2r  ) 1 20,5:0,550 

;:o  mi  = 3 

GO  TC  900 
300  Ml=3 

IF(NFSWl)  22,  25,  2? 

-»?  WRITE  (6  , 11) 

310  DO  820  11  = 1,  K 

J’C  WRIT  F (6 , A ) :H(II, J J)  , J J=1,N) 

830  WR1  TF  (6 , 12)  CELTft,*,GS 
340  IF(NFSWl)  2^,25, '>'■ 

24  WRITE (fc , 10 ) <X(I) , 1=1, N) 

GO  TC  (95  0,8  -0,850  ) , Ml 
8 r 0 WRITE  (6  , 13)  (G  (I ) , 1 = 1,  N) 

350  WRITE (6,9) 

GO  TC  (120,1 20,900) , Ml 
900  IF  (N.GT.l)  C4LL  ORDER 
CALL  FCN (N,G,F,X,M1) 

30  RETURN 
3 FORMAT  (6F12 . 5) 

7 FORMAT  (1H03E14. 3) 

8 FORMAT  (4H0  IT  I*+,rH  ST**1  14,  AH  -=*14.5) 

° FORMAT  (20M0 - - 

10  FORMAT ( 3H  OK  = 3 E14 • * / ( 3H0  6E1A.5)) 

11  FORMAT  ( 13M0F  INAL  V 4LUE*/13M0  ERRD  ? M1TRIX) 

12  FORMAT  (740DFLTA=E14. 5,-  H F=E1h.5,5M  GS  = E1^.5) 

13  FORMAT (3H05=oE14. 3/( 3H0  8E1A.5)) 

END 


SUBROUTINE  READY 

COMMON  /AGO/  H(20,  71)  ,M20)  ,G(21>  , 5 ( 2 0 ) , < 3 ( 2 0 ) , C,P  < 20  ) 

1 GB(?0)»e'»Gc»c'L»5L*rP*f>SP*T0»?,'3|A»GSS»r0»GTP*FBf&TTtGSr'» 

2 OElTA,E,N,MS,IT , L, IHPLD120) ,N?S*1,  NSSW2 

?no  i = i 

CALL  MATMPY ( N,N, H, G, 5 ) 

00  205  1=1, N 
S(1)=-S(I) 

CALL  MATMPY ( 1,N,S,S,GS) 

IF  ( 0 S *E)  21  C,2’+C,  240 
210  L=  2 

EL  = 2. 0 
T 0 = EL#F/G3 

IF  (7  0 *EL ) 217,  213,  212 

212  EL=-TO 
9i3  SL  = -GS 

DO  215  1=1, N 
215  XP(I)  =X(I)  F L* S ( I ) 

CALL  FCN ( N , G°  , co , X P, ? ) 

CALL  MATMPY  ( 1,N,  S,  G3,  G^P) 

IF  C-GSP)  240,240,  *’20 
270  IF  (F-FP)  2*0 , 24(  , ?9* 

275  L = 3 

IF  ( NFS  Wl)  100  ,101,100 

100  WRITE  (6,1) 

101  FR  = F F 

DO  220  1=1, N 
GB ( I ) =GP( I ) 

T (I) =XP(I ) 

27 0 CONTINUE 

IF  (EL-2.0)  240,2*5,740 

2*5  L=** 

DELTA  = DELTA  4-  delta 
T0=1 .0/SL 
2r*  0 RETUFN 

1 FORMAT (10H3UNOERSHOT) 

END 


SUBROUTINE  MATMPY (M,N,W,G»S) 
DIMENSION  H (20, 70  ) , S ( 70)  ,$ ( 20) 
TOO  00  720  11=1,  M 
S(II)=0.0 
00  720  JJ=1  ,N 

r> o s<ii)  =h(  jj,  id *r,(  j j)  4^(11) 

T*0  RETURN 
ENO 


SUBROUTINE  £H 

COMMON  /A 00/  H(7C,  20)  ,v(?0)  ,G(?0>  ,5(20)  ,XO(20>  ,GP(2C)  *7(20)  , 

1 G«M20)  ,r,GS,EL,SL,cP,r,SP,TO,Z,‘).  A,GSS,rO,GTP,ErfGTf  ,GSO, 

2 OELTA,P,N,MS,IT,L,IHOLD(20),NSSn,  'JSSW? 

370  L*  1 

7=3.C/EL*(F-FP)«-GS+GSP 
0 = ARS  ( 7* SORT  ( 1.0-(GS/?)*(GSP/”) ) ) 

A=  (0-7 ♦GS®)  / (0*0-SS*GSP) 

T0  = EL/3.3»  ( 0 40+?  *3  S3)  ♦/'♦A 
FO=FF-TO 

CALL  MATMPY <N,N, p, G°,T) 

PO  3 C 5 1=1, N 

305  T(l)=(GSp/3L)‘S(I) - r ( T > 

CALL  MATMpY  (t  ,N,  T,  G3,  G^P) 

IF  (TO+TC+GTF)  315,310,310 
310  DO  312  1=1, N 
312  T(I)=XP(I> *A* (X(i> -XP(T>  ) 

GO  TO  3^0 

315  IF  ( F *F  »GT  ° ) 310,720,320 

7?  C 00  322  1=1, N 
32?  TCI)  =T  (I)  +XP  (I) 

CALL  FCN(N,GT,cP,r ,?) 

IF  (FO-FB)  310,325,725 
325  L=  3 

IF(NESWl) 100 ,101, 100 

100  WRITE  (6,1) 

101  PO  327  1=1, N 
377  S(I)=T(I)-XP(I) 

CALL  MATMPY ( 1,N,S, G3,GTT) 

GTT=C-TT-GTP 
IF  ( GTT)  3-0,330,330 
330  L = 2 

GSS=CTT 
SL=-GTP 
EL  = 1 . 0 
3V  0 RETURN 

1 FORMfT (9HO^ICOCHcr ) 

END 
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SUBROUTINE  FIRE 

COMMON  /A33/  H(2C, 20) , X (20) ,G(20) , >( 2 0) ,<3< 2G> ,Gc ( 20) ,T ( 20 ) , 

1 G8(20)  ,F,GS,  EL,3L,FP,GSP,T0,?,3,A,3SS,F0,GTO,FF,GTT,GSB, 

2 OELTA,E,N,MS,IT,L,  IH0LP(2G ),N'dHl,NSSW2 
FOUIVALENCE  (TEMr,r,TT) 

400  L = 1 

TEMP=A/(1.0-A) 

CALL  FCN ( N * GO  » * T » ’) 

CALL  MATMPY (1»N»~»GT»GSR) 

T 0 = F 

IF  (TO-FP)  403,  U33,  402 
un?  T0=FF 

’♦03  IF  (T0-F‘3  + E)  415,  40 F » *05 
4T5  GSS=  C *0 

T0=Gf8*  (TEMP  -1 .0/rcvjp) 

IF  (ABS  (TO)  -0)  4-^0,41^,410 
♦ lO  L = 2 

GO  TO  *40 
*1'  L=3 

IF(FP-F)  425  ,420,420 
CO  IF  (N^SWl)  100,  1 Cl,  1 00 

100  WRITE  (6,1) 

101  FL=(1.0-A)*EL 
Fp  = FB 
GSP=CSR 

00  422  1=1, N 
XP(I)sT(I) 

GP<I>*GB<Z> 

42?  CONTINUE 
GO  TC  *40 

425  IF(NSSWl) 200,201,200 
200  WRITF (6,2) 

?n 1 FL=EL*  A 
F=F8 
GS=GSB 

00  4 27  1=1,  N 
X (I) =T (I) 

G ( I) =GP ( I ) 

CONTINUE 
GO  TC  440 
'♦T0  GSS=GSS*T  0 
00  435  1=1, N 

4*5  G(I)  =(GB(I)  -G(  I)  ) * TEMP*  (GPU  ) -G2(  l ) ) /TE*° 

440  RETURN 

1 FORMAT (lCHO^OVE  LEFT) 

2 FORMAT ( 11H0MCVE  RIGHT ) 

END 
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SUBROUTINE  DRESS 

COMMON  /A 00/  HC2C,  2J) , v(20) ,G(20)  , 3(  20), *3(20), GP ( 20 ) , T < 23 ) , 

1 GB(20),F,  GS  ,EL,  SL,  r P , GS P, T P , 7 , 0 , A , 3 SS,  c0  , DTP,  ^ , f.TT  , GSB , 

2 OELTA,E,NtMS,Il ,L,IHOLO(?0) ,NS S* 1 , N SSW2 
^00  GO  TO  (505,  E20,S'»:  , 32E)  , L 

50*  CALL  MATMPY  (N,N,  M,  3,  X) 

CALL  MATMPY ( 1 ,N, X, G, T 0) 

IF  (TO-GSS**  2/SL-E)  51r,5l0,510 
310  00  512  II  = 1 » N 
DO  512  JJ=1,N 

512  H(II,JJ)  = H(II,  JJ)  -X(IH'X(JJ)/T1 
oelta=oelta*(fl*3SS/to> 

T0=EL/GSS 
GO  TO  525 

5 lr  IF(NSSWl)  20  0 ,520,20  1 
200  WRITE  (6, 1) 

5*»0  OELTA=DElTA«  (EL*  RL  /GSS) 

T 0 = EL /GSS- 1 . C/SL 
f?5  00  527  11  = 1 , N 
00  527  J J = 1 , N 

327  H(II,JJ)=H( II, JJ) ♦ TO*  S ( 1 1 ) *S (JJ) 

510  IT=IT*1 
F=FB 

IF(NSSWl) 100,  101,100 

100  WRITE  (6  ,■* ) I T,  MS  , r , GS 

10 1 00  572  1=1, N 
G(I) =GB(I) 

XU)  =T  (I) 

CONTINUE 

IF(NSSWl) 535, 5*0, 575 
5T5  WRIT  F (6,2)  (X(I),I  = 1,N) 

WRITE (6,3)  DELTA 
540  RETURN 

1 FORMAT (BH03 OLl NE* R ) 

2 FORMAT(3HOX=  ,dE14.5/(3M0  ,ttEl4.5>) 

3 FORMAT (7H0DELTA=, *14. 5/2GH0- ) 

A FORMAT (4H0 1 T ,14, 7H  STEP  ,I4,4H  ?s,El4.5,5H  5S=,E1m.5) 

END 

SUBROUTINE  0R3FR 

COMMON  /A  00/  HH( 20 , ?0 ) , PER( 2 0) , G? ( 2 ) ) ,S(?0),XP(?0),GP(20),T(20), 

1 GB (20) ,F,GE,EL,SL,rP,GSP,T0,Z,0,A,GSS,r0,GTP,Ep,GTT,GSB, 

2 OELTA,E,N,MS,IT,L, I H OLD ( ?0 ) , NS 3 , N SSW 2 
JC  = N - 1 

DO  100  1=1, JC 
K = I ♦ 1 
00  ICO  J=K, N 

IF  (PER  (I ) -P  FR  (J ) ) 10  0,100,60 

SO  SAVE  = PER(I) 

PERU)  = ®ER(  J ) 

PER(J)  = SAVE 
SAVE  = HH ( I , I ) 

HH  (I  , I ) = HH  ( J , J ) 

HH  ( J , J ) = SAVE 
ISA  V = IHOLO(I) 

IHOLO(I)  = IHOLO(J) 

IHOLO(J)  = ISA  V 
100  CONTINUE 
RETURN 
ENO 


117 


— 


Appendix  1 

Miscellaneous  Plots 
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FIG  25.  RR  OEF  * 11  (UNOXIOIZEO)  0.6  KEV/CHRNNEL  40.000  StC 


FIG  26.  Rfl  DEF  * 12  (UNOXIDIZEO)  0.6  KEV/CHflNNEL  40.000  SEC 


F 10-  27.  Rfl  OEF  * 13  IUN0X101ZE0)  0.6  KEV/CHflNNEL  40.000  SEC 


o 


,oi»  sinhoo 


FIC  30.  Rfl  DEF  » 16  (UN0XI01ZE0!  0.6  KEV/CHRNNEL  40.000  SEC 


Fig.  31.  Randomized  Gaussian  Peaks.  Separation  ® 100  Channels 
Heights  are  10000,  20000. 


Fig.  33,  Randomized  Gaussian  Peaks.  Separation  = 100  Channels 
Heights  are  1000,  2000. 


Fig.  36.  Randomized  Gaussian  Peaks.  Separation  ■ 50  Channels 
Heights  are  100,  200. 
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